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d 

O^. Ca^"^ currents in neurons and muscle cells have been clas- 

sified as being one of 5 types of which four, L, N, P/Q and 
R were said to be high threshold and one, T, was designated 
low threshold. This review focuses on quantitative aspects 

t:;J-' I of L-type currents. L-type channels are now distinguished 

■^ ' according to their structure as one of four main subtypes 

Cay 1.1 - Cavl.4. L-type calcium currents play many funda- 

LT* I mental roles in cellular dynamical processes including pace- 

making in neurons and cardiac cells, the activation of tran- 
scription factors involved in synaptic plasticity and in im- 
mune cells. The half-activation potentials of L-type currents 

'k> I (IcaL) have been ascribed values as low as -50 mV and as 

^ ' high as near mV. The inactivation of IcaL has been found 



to be both voltage (VDI) and calcium-dependent (GDI) and 
the latter component may involve calcium-induced calcium 
release. GDI is often an important aspect of dynamical mod- 
els of cell electrophysiology. We describe the basic compo- 
nents in modeling IcaL including activation and both voltage 
and calcium dependent inactivation and the two main ap- 
proaches to determining the current. We review, by means of 
tables of values from over 65 representative studies, the var- 
ious details of the dynamical properties associated with IcaL 



that have been found experimentally or employed in the last 
25 years in deterministic modeling in various nervous system 
and cardiac cells. Distributions and statistics of several pa- 
rameters related to activation and inactivation are obtained. 
There are few reliable experimental data on L-type calcium 
current kinetics for cells at physiological calcium ion con- 
centrations. Neurons are divided approximately into two 
groups with experimental half-activation potentials that are 
high, ~ -18.3 mV, or low, ~ -36.4 mV, which correspond 
closely with those for Cavl.2 and Cavl.3 channels in physi- 
ological solutions. There are very few experimental data on 
time constants of activation, those available suggesting val- 
ues around 0.5 to 1 ms. In modeling, a wide range of time 
constants has been employed. A major problem for quan- 
titative studies due to lack of experimental data has been 
the use of kinetic parameters from one cell type for others. 
Inactivation time constants for VDI have been found ex- 
perimentally with average 65 ms. Examples of calculations 
of IcaL are made for linear and constant field methods and 
the effects of GDI are illustrated for single and double pulse 
protocols and the results compared with experiment. The 
review ends with a discussion and analysis of experimental 
subtype (Cavl.l - Cavl.4) properties and their roles in nor- 
mal, including pacemaker activity, and many pathological 
states. 
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1 Introduction 

1.1 Perspective 

The passage of ions across cell membranes, and within cells, is of 
fundamental importance in determining the electrophysiological re- 
sponses of nerve and muscle cells. Such responses are manifested 
ultimately in the functioning of the nervous and muscular systems, 
including organs of crucial biological importance such as the heart 
and brain. In the late 19th and early 20th centuries, key discoveries 
were made and biophysical theories proposed concerning such ionic 



currents, for example by Nernst (1889), Planck (1890) and Bern- 
stein (1902). With new electrophysiological recording techniques, 
many advances were made in the 1940's and 1950's by, amongst 
others, in alphabetical order, Eccles, Hodgkin, Huxley and Katz - 
see Huxley (1959) for a summary. In the 1970's and 1980's, much 
additional insight was obtained when recordings were made of cur- 
rents through single ion channels, notably by Neher and Sakmann 
(Hamill et al., 1981). In the last 20 or so years there has been an 
enormous number of discoveries concerning the factors which deter- 
mine ionic current flows in neurons and muscle cells. The present 
review concerns modeling aspects of the class of calcium currents 
called L-type, which, as will be seen below, have many consequences 
beyond electrophysiology. 

For graphic but brief historical accounts of calcium current dis- 
coveries see Tsien and Barrett (2005) and Dolphin (2006). Ac- 
cording to the former review, "...it is apparent that Ca^+ channels 
have reached the forefront of the field of ion channel research. . . due 
to their vital role in cellular signaling, their diversity, and great 
susceptibility to modulation. . . ". Records of the first single chan- 
nel recording of currents identified as being L-type were given in 
(Nowycky et al., 1985). More recent single channel recordings are 
in Cens et al. (2006), where a comparison of results for Ca^"*" and 
Ba^"*" as charge carrier is shown, and Schroder et al. (1998), where 
the much greater magnitude of L-type currents in failing heart are 
compared with those in normal heart. 

The principal motivation for the analysis and quantitative mod- 
eling of L-type calcium currents is that they occur in most nerve and 
muscle cells. They often play basic roles in pacemaker activity (see 
Section 5.2) and more generally in regulating spike frequency by in- 
ducing afterhyperpolarization, as for example in the hippocampus 
by coupling to SK channels (Tanabe et al.,1998). Wu et al. (2008) 
showed that L-type Ca^^ current in CAl pyramidal cells, by cou- 
pling to delayed rectifier potassium channels (K^7.x), can give rise 
to long-lasting changes in adaptation. 

Comprehensive models of nerve cells may include spatial vari- 
ations or not, but in either case the minimum number of current 
components is at least 10 and amongst these there should or will 
usually be included several Ca^"*" currents. If they are included in a 



model, L-type currents require a careful treatment and our aim here 
is to attempt to summarize several details of their basic properties 
and modeling which have been employed for many kinds of nerve 
and muscle cell. 

1.2 Ion channels and neurons 

Many protein molecules are embedded in the cell membranes of 
neurons. Some of these molecules are receptors for the main central 
neurotransmitters glutamate (excitatory) and GABA (inhibitory), 
as well as transmitter/modulators such as noradrenaline, dopamine, 
and serotonin all of which are released from vesicles in response to 
signals arriving at synapses where neighboring cells make close con- 
tact. (See for example Cooper et al. (2003) for an introduction to 
basic neurochemistry and neuropharmacology.) Of particular im- 
portance in determining the way in which a neuron behaves in re- 
sponse to electrical and chemical stimuli are other protein molecules 
which serve as entrance and exit pathways for electrically charged 
ions. Such molecules are called ion channels. 

If an ion channel is relatively more selective for a certain kind 
of ion, for example, sodium, Na"*", then it is called a sodium chan- 
nel. The most commonly occurring cation channels in neurons are 
sodium, potassium and calcium. Such ion channels may be open or 
closed, which means they may or may not permit the passage of ions 
through them. One of the chief consequences of the ionic currents 
which flow through such channels when they are open is the alter- 
ation of the electrical potential difference across the cell membrane. 
In the resting state, this potential difference is in the approximate 
range from -80 mV to -50 mV. An inward flux of positive ions such 
as Na"*" or Ca^"^ leads to a diminution of the potential difference, 
called depolarization or excitation, which, if sufficiently strong, may 
give rise to an action potential or spike. When the membrane is 
depolarized, ion channels such as those of Ca^"*" undergo conforma- 
tional changes which lead to their opening for certain time intervals. 
Hence they are called voltage-gated ion channels. 



1.2.1 Activation and inactivation 

The process of opening the channel is called activation. However, 
channels may be in several different states because often there is 
also a process of inactivation which is not simply the cessation of 
activation. For a channel to be conducting, the inactivation compo- 
nent must be switched off and the activation component switched 
on. If we denote the probability of activation being on by m (for 
example) and the probability of inactivation being off as h, then Ta- 
ble 1 gives the various channel states and their probabilities. It will 
be seen that the sum of the probabilities in the third column of the 
table add to unity, as they must. For more details, see, for example, 
texts such as Levitan and Kaczmarek (1997) or Koch (1999). 



Table 1: Channel configurations with activation and inactivation 



Activation process Inactivation process Probability 



Conducting state 



Off (deactivated) 
On (activated) 
On (activated) 
Off (deactivated) 



On (inactivated) 
On (inactivated) 
Off (de-inactivated) 
Off (de-inactivated) 



(1 — m)(l — h) Non-conducting 

TTlyL — llj Non-conducting 

TTlh Conducting 

(i — TTljh Non-conducting 



In the case of voltage-gated ion channels, the activation and in- 
activation probabilities usually depend on voltage and time. Hence 
if V denotes membrane potential and t denotes time, then m = 
m{V,t) and h = h{V,t). As in the pioneering work of Hodgkin 
and Huxley (1952) on squid axon, the activation and inactivation 
variables satisfy approximately differential equations whose solution 
enables one to predict reasonably accurately the ion currents flow- 
ing across nerve or other membrane. This is further elaborated on 
in Section 2.3. 



1.3 Voltage- gated calcium channels 

Calcium currents, which are found in all excitable cells, were divided 
into the two main groups of low-threshold or low-voltage activated 



(LVA) and high-threshold or high- voltage activated (HVA). The for- 
mer group contains the T-type (T for transient) and the latter group 
consists of the types L, N, P/Q and R (L for so called long-lasting, 
N, either for neither T nor L, or neuronal, P for Purkinje, and R for 
resistant). For interesting reviews of the history of the discoveries of 
these various types and the experiments that preceded them, back 
to 1953, see Tsien and Barrett (2005) and Dolphin (2006). 

Although L-type currents, the main topic of the present ar- 
ticle, were originally designated as belonging to the HVA group, 
their properties are diverse (Avery and Johnston, 1996; Lipscombe, 
2002), varying greatly with cell-type and ionic environment - see the 
data below and in particular the review of Lipscombe et al. (2004). 
They have been implicated as playing an important role in pace- 
maker activity in some neurons and cardiac cells (see for example, 
Kamp and Hell, 2000; Mangoni et al., 2003; Brown and Piggins, 
2007; Marcantoni et al., 2007; Putzier et al., 2009a, b; Vandael et 
al., 2010) and being involved in the amplification of certain synaptic 
inputs (Bui et al., 2006). 

Calcium channel molecules have up to four subunits, ai, a2-S, 
P and 7, which may exist in different forms, and which modulate 
the conductance and dynamical properties of the channel (Dolphin, 
2006, 2009; Davies et al., 2010). The ten forms of the conducting 
pore subunit, ai, lead to an expansion of the above groups (Catterall 
et al., 2005; Dolphin, 2009) to 10 main subtypes. According to 
the accepted nomenclature, L-type channels consist of the subtypes 
Cayl.l — Cavl.4. The remaining "high-threshold" currents, P/Q, N 
and R are respectively Cav2.1 -Cav2.3 and the T-current subtypes 
are Cav3.1 — Cav3.3. Of the L-channel subtypes, sub-type Ca^Ll 
is mainly found in skeletal muscle and the subtype Cavl.4 is found 
mainly in retinal cells (Lipscombe et al., 2004; Baumann et al., 2004; 
Catterall et al. 2005; Lacinova, 2005; Dolphin, 2009; Weiergraber 
et al., 2010). In cardiac myocytes and most central nervous system 
cells the subtypes are mainly Cavl.2 or Cavl.3. See subsection 5.1 
for further discussion. 

In Figure 1 is shown a recent schematic representation of the 
molecular structure of a voltage-gated calcium channel, showing the 
subunits ai, [5 and 0:2 — 5. Note that until recently the 5 subunit 
had been thought to straddle the membrane rather than project into 




Figure 1: Recent schematic of a voltage-gated calcium channel, 
taken with permission from Davies et al. (2010). The main pore 
containing subunit which straddles the membrane is ai which has 
4 juxtaposed components. The (3 subunit is intracellular and the 
a2 — S subunit projects into the extracellular compartment. In this 
example there is no 7 subunit. 

the extracellular compartment (Bauer et al., 2010; Dolphin, 2010). 
Figure 2 schematically depicts some of the structures involved in the 
important phenomenon of calcium-dependent inactivation, which is 
elaborated upon in Section 2. 

Since there are several different types of calcium currrent, and 
in particular high threshold ones, whose current-voltage relations 
are not completely distinguishing, in order to identify them, phar- 
macological methods are used. Each channel type/subtype, acts as 
a receptor for specific molecules, most often toxins, leading to re- 
duced calcium current. In some cells, R-type currents, once thought 
not to have a specific blocker and which are borderline between 
HVA and LVA, have been more recently found to be blocked by 
the tarantula toxin SNX-482 (Li et al.,1997; Newcomb et al.,1998). 
Table 2, adapted from Striessnig and Koschak (2008) and Catter- 
all et al. (2005), summarizes the main agents used for blocking or 
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Figure 2: Sketch representation of an L-type Ca^^ channel depict- 
ing the mechanism of GDI. Calcium ions are depicted as filled red 
circles. The calcium binding protein calmodulin is shown tethered 
to one of the intracellular motifs alongside the EF hand. For the 
role of calmodulin in inactivation, see Section 2.1. This figure has 
been composed from figures in Bodi et al. (2005), Dolphin (2009), 
Davies et al. (2010) and Doan (2010). 

partially blocking the various types/subtypes. Note, however, that 
dihydropyridines, phenylalkylamines and benzothiazepines block all 
L-type channels to some degree. 



It is noteworthy that one of the most important areas in which 
the quantitative study of L-type (and other) Ca^+ currents plays 
a fundamental role is in the electrophysiology of the heart. In 
medicine, L-type calcium currents are the prime target of impor- 
tant calcium channel blockers such the benzothiazapine diltiazem, 
the phenylalkylamine verapamil and the dihydropyridines such as 
nifedipine, which are used to treat, inter alia, hypertension, angina 
and cardiac arrhythmias. 

Within the above subtypes, various configurations of other sub- 
units lead to channels, with quite different properties (Dolphin, 
2009; Catterall, 2010). Thus, L-type Cavl.3 variants may have 
differing current magnitudes and pharmacology (Andrade et al., 
2009). One cannot therefore ascribe definite parameters in the dy- 
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Table 2: Summary of pharmacological agents used 
to identify various types/subtypes 



Type 


Subtype by 
ai subunit 


Agent (Example) 


L 


Cavl.l 


Dihydropyridines (Nifedipine) 




Cavl.2 


Phenylalkylamines (Verapamil) 




Cavl.3 


Benzothiazepines (Diltiazem) 




Caa-4 




P/Q 


Cav2.1 


w-Agatoxin IVA 


N 


Cav2.2 


w-Conotoxin-GVIA 


R 


Cav2.3 


SNX-482 


T 


CavS.l 
Ca^3.2 
Cav3.3 


Kurtoxin 



namical description of L-type calcium currents based on the sub- 
types Cavl.l — Cavl.4. For example, the (3 subunit regulates cur- 
rent density (defined, for example, as current per unit area) and 
the properties of activation and inactivation. There are four forms 
of the (3 subunit, [3^ and /34 being highly expressed in the brain 
(Dolphin, 2003). 

1.3.1 A brief background and useful references 

It has become apparent that quantitative descriptions of L-type cal- 
cium currents are often a fundamental component in computational 
modeling of neuronal and muscle cell dynamics. The first descrip- 
tions of L-type current dynamics in such a context appear to be 
those of Belluzzi and Sacchi (1991) and McCormick and Huguenard 
(1992), although there had been many models containing calcium 
currents, prior to the division into the above types, for example 
for cardiac cells. In fact Belluzzi and Sacchi's (1991) model for a 
sympathetic neuron in the rat superior cervical ganglion described 
a high-threshold current which had an inactivation similar to the 
L-type of Fox et al. (1987) but an activation similar to the N-type 
current. 

For an introduction to fundamentals about ion channels see, 
for example, Levitan and Kaczmarek (1997). Additional reviews 

11 



on voltage-dependent Ca^^-channels pertaining to their molecular 
structure, nomenclature and function, including modulation and 
regulation are contained in, for example, in chronological order, 
Tsien et al. (1988), Bean (1989), Tsien and Tsien (1990), Anwyl 
(1991), Catterall (1995), De Waard et al. (1996), Jones (1998), Hof- 
mann et al. (1999), Catterall (2000), Ertel et al. (2000), Kamp and 
Hell (2000), Lacinova (2005), Zamponi (2005) and Catterall (2010). 
Reviews addressing specific related topics include those on synaptic 
transmission (Meir et al., 1999; Neher and Sakaba, 2008), T-type 
currents (Perez- Reyes, 2003; Cueni et al., 2009), calcium-dependent 
inactivation in neurons (Budde et al., 2002), dynamics of calcium 
signaling in neurons (Augustine et al., 2003), models of calcium 
sparks and waves (Coombes et al, 2004), L-type currents in the 
heart (Bodi et al., 2005), the role of calcium currents in circadian 
rhythms (Brown and Piggins, 2007), calcium release and the roles 
of ryanodine receptors in heart and skeletal muscle diseases (Zalk, 
2007), Ca^"^ channels in chromaffin cells (Marcantoni et al., 2008) 
and calcium dynamics in relation to absence epilepsy (Weiergraber 
et al. 2010). Calcium ion infiux through L-type channels leads via 
various signaling pathways to the activation of transcription factors 
such as CREB and hence the expression of genes that are essen- 
tial for synaptic plasticity and other important cellular processes 
(Dolmetsch et al., 2001; Hardingham et al., 2001; Mori et al, 2004; 
Power and Sah, 2005; Satin et al., 2011). The role of calcium ion in- 
fiux through L-type channels in immune cells has also been recently 
reviewed (Suzuki et al. 2010). It was also recently demonstrated 
that calcium currents of L-type, together with R-type, are involved 
in the activation of SK channels, which attenuate excitatory synap- 
tic transmission in pyramidal cells of the medial prefrontal cortex 
(Faber, 2010). 

Mathematical modeling of the electrophysiology of cardiac cells 
has a history spanning nearly the last 50 years - see for exam- 
ple DiFrancesco and Noble (1985) - and reviews by Noble (1995), 
Wilders (2007), Fink et al. (2011) and Wilhams et al. (2010). 
Brette et al. (2006) reviews background physiology and biophysics 
of calcium currents in cardiac cells. The modeling has taken into ac- 
count details of the dynamics of several ionic currents, particularly 
calcium currents (as for example in Luo and Rudy, 1994). For a 
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discussion of the dynamics underlying the large difference between 
the relatively short action potential duration in neurons and heart 
cells, see Boyett et al. (1997). In ventricular cardiomyocytes L-type 
calcium currents play a pivotal role (Benitah et al., 2010) although 
T-type currents, important in development and in some patholo- 
gies (Ono and lijima, 2010) are also present (Bean, 1989) and play 
a role in pacemaker activity (Bers, 2008). L-type channels, located 
in the sarcolemma, are involved in arrythmias, and are activated by 
depolarization and inactivated by both voltage-dependent and in- 
tracellular Ca^^ mechanisms. The inactivation of L-channels, which 
is elaborated upon below, is complex, not only in cardiac cells, and 
not completely understood (Imredy and Yue, 1994; Findlay et al. 
2008; Grandi et al., 2010). See also Lacinova and Hofmann (2005) 
and the review of mechanisms of calcium and voltage-dependent in- 
activation in Cens et al. (2006). The website of David Yue (Calcium 



signals lab) http://webl.johnshopkins.edu/csl/ contains a wealth of 



information about calcium dynamics. 

In the electrophysiology of neurons and muscle cells, calcium 
currents usually have crucial roles. Hence there is the need to rep- 
resent as accurately as possible, in mathematical or computational 
modeling, calcium entry into cells along with related processes, such 
as buffering and pumping. The latter two topics are not explored 
in depth here, but see for example Standen and Stanfield (1982), 
Blaustein (1988), Tsien and Tsien (1990), McCormick and Hugue- 
nard (1992), Tank et al. (1995), Lindblad et al. (1996), Blaustein 
and Lederer (1999), Stokes and Green (2003), Shiferaw et al. (2003), 
Rhodes and Llinas (2005), Roussel et al. (2006), Higgins et al. 
(2007), Friel and Chiel (2008) and Brasen et al. (2010). 

2 Quantitative description of L-type Ca^+ 
currents 

We are here concerned with deterministic approaches that have been 
employed in the quantitative description of L-type Ca^"*" currents. 
Such descriptions are probably adequate in many preparations al- 
though it has not been found to be the case for high rates of action 
potentials in cardiac myocytes, where a many-state Markov chain 
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model, used in conjunction with the constant field method (Maha- 
jan et al, 2008), was found to be more accurate. See also Jafri et 
al. (1998), Sun et al. (2000), Fink et al. (2011), Grandi et al. 
(2010) and Williams et al. (2010) for more details on the Markov 
chain approach for determining L-type channel-open probabilities. 
Destexhe and Sejnowski (2001) contains a brief introduction to such 
models, especially for voltage-dependent sodium channels. As far 
as can be discerned, until now there have been no Markov models 
of L-type calcium channels in non-cardiac cells. 

2.1 General description: VDI and GDI 

Whereas activation kinetics of L-type channels are voltage-dependent 
(Budde et al., 2002; Lacinova, 2005), inactivation has often been 
found to have dependence on both voltage and the internal "con- 
centration" of calcium ions, as discovered by Brehm and Eckert 
(1978) in Paramecium and Tillotson (1979) in molluscan neurons. 
The two components are called GDI (calcium-dependent inactiva- 
tion) and VDI (voltage-dependent inactivation). A hallmark of GDI 
is that it has a maximal effect where the relevant calcium currents 
are themselves maximal. Other markers are the reduction in GDI 
if Ba^+ is the charge carrier or if the buffers EGTA or BAPTA are 
present in the pipette (see Budde et al., 2002, for a discussion). 

The relative contributions of GDI and VDI vary widely amongst 
cell types and modeling these contributions is not always as straight- 
forward as in the simplified scheme presented below. Sometimes 
GDI may be absent, especially in the case of Gavl.4 channels (Koschak 
et al., 2003; Baumann et al. 2004, Singh et al., 2006; Wahl-Scott 
et al., 2006; Liu et al, 2010). Indeed, often GDI is not included 
in modeling and often there is no inactivation at all, neither volt- 
age nor calcium dependent, but it is not clear if this renders some 
calculations inaccurate. 

In Tables Al.l-Al.4 of the Appendix, there are given data for 
about 65 studies, which involved L-type (or undefined high thresh- 
old calcium) currents. This survey is not exhaustive but represen- 
tative of works from 1987 to the present. About half of these are 
experimental and about half concern mathematical modeling. Some 
salient aspects of the experimental studies are discussed in Section 
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3. Here we note that for the 24 neuronal modehng studies, 16 in- 
cluded no inactivation, 3 included VDI only, 3 included GDI only 
and a further 2 included both VDI and GDI. It is not clear if in- 
activation was omitted in several cases simply because the models 
employed data extrapolated from other cell types or the inactiva- 
tion was considered to be unimportant or extremely slow. For the 7 
models for cardiac cell L-type calcium currents, 2 studies included 
VDI only and 5 included both GDI and VDI. It is claimed (De 
Waard et al., 1996) that inactivation of L-type currents is generally 
slow in neurons and secretory cells but more rapid in cardiac cells. 
The data in Table A3 indicate that this may not always be the case, 
but the definition of the term rapid is no doubt flexible. 

Generally, GDI has been cited as being more rapid than VDI 
(Budde et al., 2002; Lacinova and Hoffman, 2005; Grandi et al., 
2010). The molecular basis of the dynamics of GDI has been an 
active and fascinating research area in the last two decades (Imredy 
and Yue, 1992,1994; Lee et al., 1999; Peterson et al., 1999; Qin 
et al., 1999; Ziihlke et al., 1999; Erickson et al., 2003; Soldatov, 
2003; Bazzazi et al., 2010). It was found that calmodulin, which is 
tethered to the channel, must bind Ga^"*", as depicted in Figure 2, 
whereupon a configurational change takes place resulting in inac- 
tivation. Grump et al. (2011) have reported that calmodulin and 
Ga^"*" can compete to limit GDI in Gavl.2. The roles of the /3 and 
a25 subunits in determining the properties of Gavl.2 channels have 
been discussed in Ravindran et al. (2008) and Ravindran et al. 
(2009). 

There have been developed, since the original discovery of GDI, 
two main modeling ideas concerning the spatial distribution of cal- 
cium ions which participate in the inactivation process. One, called 
shell theory and posited originally by Standen and Stanfield (1982), 
is that there is a region of some 100 nm depth inside the cell mem- 
brane where calcium ions may accumulate, giving a concentration 
Ga* much higher than in the remaining cytoplasm and that it is Ga* 
which should be used to determine the rate of GDI. This approach 
is widely used in neuronal modeling, if indeed GDI is included. The 
second approach, called domain theory (Sherman et al., 1990), takes 
account of the calcium ions just inside the pore where they have en- 
tered the cell. In fact, according to Imredy and Yue (1994) and 
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Cens et al. (2006), inactivation (GDI) can be induced in a sec- 
tion of membrane which contains a single L-type channel. However, 
in some cells, in a region with many channels, averaging over all 
channels should make the often used and simpler shell approach a 
reasonable approximation. 

Special consideration is made for some cardiac and other mus- 
cle cells where calcium entry, principally through L-type channels, 
leads to Ca^^-induced calcium release (CICR) of stored calcium via 
ryanodine receptors and results in rapid GDI. Many articles have 
addressed the modeling of CICR and the geometrical details and the 
formation of localized increases in Cai, called calcium sparks, and 
their triggering of a large local increase in Ca^ (see for example and 
references therein, Shiferaw et al., 2003; Soeller and Cannell, 2004; 
Koh et al., 2006; Bers, 2008; and Groff and Smith, 2008). Many of 
these approaches employ Markov models (Hinch et al., 2004; Shan- 
non et al., 2004; Greenstein et al. 2006; Faber et al., 2007). Scriven 
et al. (2010) have provided detailed information on the geometry 
and numbers of Cavl.2 channels and ryanodine receptor clusters in 
rat VM which are useful in modeling Ca^''' dynamics. 

Ryanodine receptors may also be coupled to L-type calcium 
channels in some neurons so that CICR plays a role. A comprehen- 
sive review containing many examples of CICR in neurons, including 
that involving L-type currents in the hippocampus, was compiled by 
Verkhratsky (2005). Coulon et al. (2009) report the occurrence of 
both low threshold and high threshold calcium currents in connec- 
tion with CICR in the thalamus. In bullfrog sympathetic neurons, 
Albrecht et al. (2001) showed that small elevations of Cai evoked 
by weak depolarization lead to Ca; accumulation by the endoplas- 
mic reticulum, and that Caj accumulation became stronger after 
inhibiting CICR with ryanodine. A mathematical model was pre- 
sented in support of these results. In another interesting related 
study, Hoesch et al. (2001) concluded that caffeine is a reliable 
agonist for CICR in rabbit vagal sensory neurons, but that caffeine- 
activated rises in Caj in nerve cells could not be attributed solely to 
release from intracellular stores. See Tsien and Tsien (1990), Friel 
and Tsien (1992), Chavis et al. (1996) and Ouyang et al. (2005) 
for further examples and discussions of neuronal CICR. 
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2.2 The basic model for L-type Ca^+ current 

In the following, the membrane potential is V, the internal calcium 
concentration is Caj, the external calcium concentration is Cao and 
t is time. All the deterministic formulations of the L-type calcium 
current employed in modeling to date are included in the general 
form 

IcaL = m^i(\/,t)/iP^(\/,t)/(Cai,t)F(\/,Cai,Cao), (1) 

where m{V, t) is the voltage-dependent activation variable, h{y, t) 
is the voltage-dependent inactivation variable and /(Cai,t) is the 
(internal) calcium-dependent inactivation variable. The factor F 
contains membrane biophysical parameters and, as described be- 
low, is of the Ohmic (or linear) form (as in (15)) used in the origi- 
nal Ho dgkin- Huxley model, or the constant-field form (as in (17) or 
(19)), often called the Goldman-Hodgkin-Katz form. The values of 
Pi and p2 are ideally dictated by best fits of current- voltage relations 
to experimental data. The power pi to which m is raised, is about 
equally frequently pi = 1 or pi = 2, with invariably pi = 1 for car- 
diac cells. For skeletal muscle, before the L-type was distinguished, 
the value pi = 3 was employed (Standen and Stanfield, 1982). De- 
tails are given in Tables Al.l- A1.4 in the Appendix. The value of 
P2, if indeed VDI is included, is invariably p2 = 1. (Putting P2 = 
implies no VDI.) The notation fca is often used for / and / often 
used for h, but the present notation avoids excessive subscripts in 
subsequent formulas. If other calcium currents were under consider- 
ation, the variables for L-type current might be usefully written as 
m^, hi and fi. In some reports, L-type Ca^"^ channels are described 
as being also permeable to Na^ and K+, the relative permeabilities 
being given for ventricular myocytes as 2800:3.5:1 (Luo and Rudy, 
1994; Faber et al., 2007) and 3600:18:1 (Shannon et al., 2004). The 
contributions from Na"*" and K"*" are evidently sufficient to bring the 
reversal potential for IcaL from its value around the Nernst poten- 
tial for Ca^"*" (about 120-150 mV), which it would be close to if 
it was through a purely Ca^"^-conducting channel, to the observed 
values of about 70 mV. However, in what follows there is no focus 
on the Na+ and K+ components of the L-type current. In model- 
ing neuronal or other cell-type dynamics, there are of course many 
other components, one of which will be the intracellular calcium 
concentration Cai(t) whose value will directly inffuence IcaL, Cao(t) 
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usually being regarded as constant. 

Each of the variables m, h and / takes values between and 1, 
inclusively, and the product m^^h^^f is interpretable as the proba- 
bility that the channel is open or conducting or equivalently gives 
the expected fraction of such channels in a large sample. As can 
also be seen from Tables Al.l-Al.4, there have been many forms 
other than what might be called the full description as in Equ. (1) 
in which there there is both time-varying voltage-dependent inac- 
tivation and time-varying calcium-dependent inactivation. Thus, 
as mentioned in Section 2.1, in some computations, there is no in- 
activation whatsoever, or there may be just one or the other of 
voltage-dependent or calcium-dependent inactivation, either with 
or without time dependence. Sometimes, if a time constant for an 
activation or inactivation variable is very small, the steady state 
value of the variable is assumed to hold instantaneously with no 
time dependence. With a few exceptions, this has been the case in 
modeling calcium-dependent inactivation. 

2.3 Activation and inactivation variables 

Most of the material in this subsection is standard but is repeated 
here for notational purposes and to make the account self-contained. 
According to the above basic model there are three dynamical vari- 
ables determining the magnitude of the current. These are the acti- 
vation variable, always assumed to be purely volt age- dependent, the 
voltage-dependent inactivation variable, and the calcium-dependent 
inactivation variable. 

2.3.1 Activation variable m and voltage-dependent inac- 
tivation variable h 

The activation variable m{V,t) is usually assumed to satisfy the 
differential equation 

dm _ rrioo - m 

It ~ 7 ' ^ ' 

where niooiV) is the steady state value m(V, oo) and r^ is a time 
constant which may depend on V. Similarly for the voltage-dependent 
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inactivation variable 

dh hno — h 

— = — . 3 

dt Th ^ ^ 

In the classical approach, the voltage dependence of the steady 
state activation variable is written in the Boltzmann form 

1 + e '5 

where V is the membrane potential. The half-activation potential is 
given by moo{V^ i ) = 0.5, and if the derivative of moo is denoted by 
m'^ then the slope factor is km = ^^, }y — y. Similarly, if a voltage- 

dependent inactivation is included then its steady state form may 
be written as 

1 + e "'^ 
where hoo{Vf^ i ) = 0.5 and the slope factor is kh = —jri—k? — r- 

'2 oo^^h,^' 

In many instances the Boltzmann forms are not used but rather 
forward and backward rate coefficients cxiV) and PiV) respectively, 
are introduced, so that 

dui 

— = am{y)(l - m) - I3m{y)m (6) 

dt 

with a similar equation for h. It is usually the case that moo(^) = 
— /f.T,\a nr\ and TrniV) = — ., ,^ -*- „ ,,,-, aud similarly for the inacti- 

am{V)+l3m(V) "'^ ' Oim{V)+l3m(V) -^ 

vation variable, but as will be seen from Table A2, variations have 
been employed. 

The Boltzmann forms for rrioo and hoo have the advantage of 
making it immediately transparent at which voltages these functions 
have their half-maximal values and over what ranges of V they are 
significantly different from zero or one. If a and /3 are given, and 
^oo(^) = (Y)j^n ly) ) then one may determine an approximating 
Boltzmann curve by graphical inspection or calculation to ascertain 
the Vi and /c-values. 

2 

2.3.2 Calcium-dependent inactivation variable /(Cai,t) 

It will be seen in the summary of various models presented below 
that the inactivation variable /(Cai,t) has entered the modeling in 
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many different forms. The most general, which contains all those 
employed, can be written as 

where /oo(Cai) is the steady state /(Cai,oo) and r/ is a time con- 
stant which may depend on Cai. Here /oo and tj are defined through 

1 

Joo ■ 



l+(t)" 
"' ^ "/[l+(t)"] *" 

where a/, with units time"^, is a constant and Kf in mM is the value 
of Cai at which the steady state inactivation has half-maximal value. 
The index n has been given the values 1 (Standen and Stanfield, 
1982), 2 (Luo and Rudy, 1994) and 3 (Fox et al., 2002). It can 
be argued that n, which is often referred to as a Hill coefficient, is 
the approximate number of calcium ions which bind to a channel 
or associated molecule to give inactivation. However this reasoning 
is open to considerable doubt (Weiss, 1997). Hofer et al. (1997) 
reported that the Hill coefficient for binding of Ca^"*" to the site 
mediating GDI was close to 1 with an "inhibition constant " of 
4/xM. In the original derivation of Standen and Stanfield (1982) it 
was assumed that the underlying reaction was in fact simply 

Ca + R^CaR, (10) 

Of 

with Kf = af/f3f. Using standard reaction rate theory (e.g. Mc- 
Conigle and Molinoff, 1989) this scheme gives the ratio of inacti- 
vated to activated as 

_Ca^ (11) 

^ Cai + Kf 

(Recall that / is the interpretable as the probability that a channel 
is not inactivated by a Ca^''"-dependent mechanism.) The alterna- 
tive formulation of the time-dependent behaviour in the case of n 
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calcium ions binding to give an inactivated channel is through the 
differential equation 

| = a;(l-/)-/3;Car/ (12) 

where Pj = af/K^ which is consistent with (7)- (9). However, re- 
finements of this basic model would be needed to incorporate other 
findings such as the dependence of GDI on voltage in cardiac my- 
ocytes (Faber et al., 2007). 

2.4 Determining the current 

Although they are well-known, we give here for completeness the 
two most generally used methods for calculating membrane ion cur- 
rents (see e.g., Tuckwell, 1988; Koch, 1999). These are the linear 
and constant field methods, of which the former is simpler. In the 22 
works on non-cardiac cells summarized in Table Al.l, which state 
the method of determining the current, 15 use the linear method 
and 7 use the constant field method. All cardiac cell studies in Ta- 
ble Al.4 use the constant field method. Although the constant field 
method is considered more appropriate, if the voltage doesn't spend 
much time above then the linear method may be sufficiently accu- 
rate, especially when there is voltage-dependent inactivation, which 
is supported by the calculations in Section 4.1. 

2.4.1 Linear method 

This method, employed by Hodgkin and Huxley (1952) in their fun- 
damental modeling of the action potential in squid giant axon, con- 
sists of multiplying the membrane conductance Gl by the driving 
potential V — Vrev,L, where Vrev,L is a reversal potential, to give 

IcaL.Lin = Gl{V — K-cd.l)- (13) 

This method is called linear (or Ohmic), as the term V — Vrev,L is 
linear in V, but not of course the conductance. The conductance 
Gl is the product of a maximal value Gi^max and one to three of 
the factors m(y, ty^, h(y, ty^ and /(Caj, t), all of which are dimen- 
sionless and take values between and 1 inclusively. That is, 

Gl = m{V, trh{V, ty^f {Ca.,, t)GL,ma.. (14) 
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Assuming a membrane area A sq cm, a channel density pi per 
sq cm and a single channel open conductance of gi^, the maximal 
conductance is GL,max = ApLQi- Then in terms of fundamental 
quantities 

IcaL,Li„ = ApLg^m{V,ty^h{V,ty'f{CB.ut){V - Vrev,L). (15) 

As can be seen from Tables Al. 1-1.4, pi and p2 are usually small 
non-negative integers between and 2 for L-type currents. If ^^ is 
in mS and voltages are in mV, then IcaL.Lin given by (15) is in pA. 
The reversal potential for an ion type is often taken to be its 
Nernst potential which for Ca^"*" is 

where Cao and Cai are external and internal concentrations, R, T 
and F have their usual meanings, RT/F having units of volts. Al- 
though Vca is sometimes taken as the value of Vrev,L, it is usually of 
magnitude about 120 mV, which is much higher than experimentally 
determined values of Vrev,L which are around 60-70 mV. In model- 
ing, using Vca thus gives substantially larger calcium currents than 
expected. 

2.4.2 Constant field method 

From the basic expression for the current density under the constant 
field assumption (see for example Tuckwell, 1988, Chapter 2), the 
L-type calcium current, in pA, through a membranous area of A 
square cm is given by the Goldman-Hodgkin-Katz flux equation 

r -—1 

IcaL,cF = APLmP'{V,t)h^'{V,t)f{Ca.ut)^. ^^""^ ~ ^""t "° , 

(17) 
where for convenience of expression we have defined the temperature- 
dependent voltage 

RT 

F is Faraday's constant (96500 coulombs/mole), Pl is the mem- 
brane permeability coefficient in cm/sec, and Cai and Cao are the 
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K) = ^^, (1^ 



intracellular and extracellular concentrations of Ca^'^ in mM. Note 
that the partition coefficient, here assumed the same for the intracel- 
lular and extracellular boundaries, is absorbed into the permeability 
coefficient. However, if the intracellular and extracellular partition 
coefficients are different, then additional multiplicative factors will 
appear with ion concentrations, as in, for example. Sun et al. (2000) 
and Fink et al. (2011). In terms of single channel properties, the 
current in /iA through an area of A cm^ is 



2V . 



lc.L,CF = ApLPlmP'iV,t)hP''iV,t)fiCaut)^^^ ^^""^ ^''°'' ""^ 



2V 



Vo 1 _ e vo 

(19) 
where p^ is the density of channels in cm ^ and P£ is the single 
channel permeability in cm^/sec. 

3 Data on L-type activation and inacti- 
vation 

A large number of sources of data for L-type activation and inacti- 
vation, including both VDI and GDI, have been considered, back to 
the first uses of the term "L-type" (Nowycky et al., 1985; Fox et al., 
1987). Data on the half activation potentials V^ i and slopes km are 
given in Tables Al.l-Al.4. In these and all subsequent tables the 
authorship of articles is given only by ffist author and date to mini- 
mize column width. In the ffist column of Tables Al. 1-1.4 are given 
carrier concentrations. M denotes a purely mathematical modeling 
study, the absence of M can denote either a purely experimental 
work, or experimental work in conjunction with modeling. The 
characterization of currents as L-type by various authors has been 
assumed to be correct. Voltage-dependent inactivation parameters 
Vi^ 1 and kh are also given if available. Table Al.l contains data for 
neurons and secretory cells, whereas Table Al.2 is restricted to car- 
diac cells. In some cases, marked with asterisks (*), parameters have 
been estimated approximately from data presented graphically, but 
the original articles should be consulted for details. Table Al.3 has 
data for neurons and secretory cells with both VDI and GDI and 
Table Al.4 presents these data for cardiac cells. In Table Al.l the 
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entry for Joux et al. (2001) for V^ i has two values as fits were 
to a double Boltzmann function. Note that the concentrations Cao 
may differ in the various preparations, (and occasionally be Bao) as 
described in the first column of Tables Al.1-1.4. For example, the 
value in the first entry (Fox et al., 1987) is Cao=10mM, which may 
explain why the value of V^ i seems rather high. In general, much 
of the variability in the data presented is attributable to various 
ionic carrier concentrations and differing subunit isoforms. In order 
to keep the Tables relatively simple, many such details are omitted, 
but these can be obtained from the original sources in most cases. 

3.1 Steady state activation moo{V) and inactiva- 
tion hoo{V) 

In the top part of Figure 3 the steady-state activation rrioo is plotted 
against V for the two lowest and two highest values of the half- 
activation potential found in a comprehensive, but not exhaustive, 
search of the literature, including all modeling and experimental 
results in the Tables Al.1-1.4. The two lowest are for midbrain 
dopamine neurons, with half- activation potentials of -50 mV ( Amini 
et al., 1999; Komendantov et al., 2004). 

The two highest are +5 mV for ventricular myocytes (Shiferaw 
et al., 2003) and +2 mV for chick dorsal root ganglion (Fox et al. 
1987). Also shown in the top part of the figure are examples of the 
rrioo curves found for the channel subtypes Cavl.l-Cavl.4. It can 
be seen that the curves for the subtypes cover or span those for the 
lowest and highest activation curves for cells except for the one for 
midbrain dopamine neurons. 

In the lower part of Figure 3 are shown the corresponding steady 
state V-dependent inactivation hoo curves for the two highest half- 
activation potentials. For the other two cases with l^ i=-50, there 
was only GDI. Also shown in the lower part of Figure 3 are the 
hao curves for each channel subtype. The curve for Cavl.2 seems 
more extreme than that for any cell. In some cases Boltzmann 
functions for maoiV) and hooiV) were not given but rather explicit 
expressions for the forward and backward rate functions amiV), 
/3m(y), ahiV) and f3h(y). Such formulas when available are given 
for all cell types in Table A3. The parameters of fitted Boltzmann 
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Figure 3: Steady state activation moo and inactivation h^o ver- 
sus membrane potential for L-type currents in neurons for extreme 
cases. Top panel. Activation curves rriooiV) are given for L-type 
currents in midbrain DA neurons, (blue solid curve, Amini et al., 
1999; red solid curve, Komendantov et al., 2004) chick DRG (green 
solid curve. Fox et al., 1987) and rabbit atrial myocytes (black solid 
curve, Lindblad et al., 1996). These are the cases with the two 
highest and two lowest half-activation potentials found in the lit- 
erature. The two lowest (midbrain DA) have V^ i =-50 mV, and 
the highest two are for rabbit atrial cell, with VI, i=-0.95 mV and 
chick DRG with V^ i=2 mV. Also given are examples of activation 
curves (dashed) for the four subtypes, Cavl.l (green, from Catter- 
all, 2005), Cavl.2 (blue, average from Lipscombe et al., 2004 and 
Catterall, 2005), Ca^LS (red, average from Lipscombe et al., 2004 
and Helton et al., 2005) and Cay 1.4 (black, mean of lowest and 
highest, from Catterall, 2005). Bottom panel. Steady state inacti- 
vation hoo{V) for cases shown in the top panel with the same color 
coding. For the two cells with the lowest two half- activation poten- 
tials, inactivation was given only as Ca^^-dependent - see Tables 
A1.1-A1.4. 
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functions for the corresponding 171^0 and h^o were given in Tables 
Al.l-Al.4. When exphcit expressions were given for the activation 
time constants Tm{V), these are given in Table A4. 

3.1.1 Distribution of half-activation potentials 

It can be seen from Tables Al.l-Al.4 that there is great variability 
in the parameters for the activation and inactivation properties of 
L-type currents in various cells. The underlying reasons for variabil- 
ity include, (a) the occurrence of different subtypes of L-channels, 
some examples of whose properties are listed in Table A6; (b) the 
various locations of the channels over the soma-dendritic surface; 
(c) the interaction of L-type channels with other neighboring ionic 
channels; and (d) the ionic compositions of the intracellular and ex- 
tracellular compartments of the cell under consideration. One other 
source of variability is that there are different methods, some more 
accurate than others, for calculating Boltzmann curves. 

The majority of data in tables Al.l-Al.4 concerns modeling 
studies and it is useful to separate out the experimental values, 
and to subdivide these into those for neurons and cardiac cells. In 
particular, the experimental results considered here are only those 
for calcium concentrations in the physiological range because high 
concentrations of either barium or calcium ions can shift the half- 
activation potential in the depolarizing direction by 10 or more mV 
(Jaffe et al., 1994; Muinuddin et al., 2004). With these restrictions 
there are 13 half- activation potentials given for neurons and 5 for 
cardiac cells. Histograms for the l^ 1 values are shown in Figure 
4. For neurons, the data are divided into a high (> —30 mV) and 
a low group, the topmost panel being for the low group. The mean 
for the neuron high group is < V^i >= —18.3 mV whereas that 
for the low group is < V^ 1 >= —36.4 mV. The mean slope factors 
for these two groups of neurons are < km >= 8.4 mV for the high 
group and < km >= 5.8 mV for the low group. These average slope 
factors are obtained on omitting the outliers 2.45 in the high group 
and 20.0 in the low group. The mean value of Vm i for the exper- 
iments with cardiac cells is -13.6 mV with a standard deviation of 
7.8 mV, and the corresponding average km is 6.84 mV with standard 
deviation 1.61 mV. It is then of interest to note that the average 
values of the half activation potentials for Cay 1.2 and Cavl.3, using 
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Figure 4: Histograms of the half-activation membrane potentials 
V^ 1 for L-type Ca^^ currents from the data of Tables Al.l-Al.4. 
The data, which are based on selected experiments (see text) are 
divided into 3 groups. Top, low neuron group. Middle, high neuron 
group. Bottom, cardiac cells. Note that the mean values for Cavl.2 
and Cayl.S, using Table A6 values with physiological calcium ion 
concentrations, are -18.7 mV and -38.9 mV, respectively. 

Table A6 values with physiological calcium ion concentrations, are 
-18.7 mV and -38.9 mV, respectively. These values correspond well 
with the mean values of the high and low groups of neurons. The 
values for cardiac cells tend to be higher than those in even the 
high neuron group, but the sample sizes are too small to draw any 
statistically significant conclusions. 

For all the 25 neuronal modeling studies the mean value of VL i 

'"'2 

is -21.0 mV with a standard deviation of 12.21 mV, and for k^ the 
mean is 7.26 with a standard deviation of 2.33. The mean modeling 
study parameters correspond roughly to the experimental values for 
the high neuron group. For the 7 cardiac cell (all kinds) modeling 
studies the mean half-activation potential V^ i is -8.8 mV with a 
standard deviation of 9.1 mV and the corresponding mean for km 
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is 6.4 mV with a standard deviation of 0.52. 

For the parameters of the Boltzmann function for the steady 
state inactivation, V^ i and kh, there are, unfortunately httle ap- 
propriate experimental data, and none which could be considered 
reliable, though many modeling studies have employed estimates 
(see Tables Al.l-Al.4) and see below. In order to provide crude 
estimates based on experiment, one may use the results of Koschak 
et al. (2001) as given in Table A6 for the differences in voltage 
between V^ i and Vf^ i for Cavl.2 and Cay 1.3 in high barium con- 
centration. This procedure yields values which are summarized in 
Table 3 for the activation and inactivation parameters. 



Table 3: Parameters for L-type steady state 
activation and inactivation in neurons 





High neuron 


Low neuron 


All models 




group (expt) 


group (ex 


:pt) 




Ka 


-18.7 


-38.9 




-21.0 


rCfYi 


8.4 


5.8 




7.26 


Vh.i 


-42.0* 


-61.9* 




-38.95 


h 


13.8* 


6.6* 




6.82 



* Estimated^- see text. 

For cardiac cells, based on the few experimental studies available, 
the mean Vu i is -40.0 mV with standard deviation 12.8, the corre- 

"'2 

sponding mean value of kh being 6.85 mV with a standard deviation 
of 2.04 mV. 

Of the 25 neuronal modeling studies in the tables, only 6 in- 
cluded VDI and these yielded a mean value of -38.95 mV for Vf^ i 
with standard deviation 14.29. For the corresponding values of the 
slope factor kh, the mean is 6.82 with a standard deviation of 2.16, 
where one value is omitted as an extreme outlier. For the cardiac 
modeling, the mean value of Vf^ i is -30.1 mV with a standard devi- 
ation of 10.9 mV and the corresponding mean for kh is 6.1 mV with 
a standard deviation of 2.00 mV. 
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3.2 Time constants rjn{V) for activation 

Some values of Tm{V) for activation, with corresponding inactiva- 
tion, if available, are reported in Table A3. This table has been 
divided into a part with experimental results and a part on values 
of which have been employed in modeling studies. As seen in Table 
A3-Experimental, there are not many explicit experimental results 
available for the time constants of activation (or inactivation) of L- 
type Ca^"*" currents. In neurons values of Tm{V) range from about 
0.5 ms in sympathetic neurons (Belluzzi and Sacchi, 1991) to 2.3 
ms in hippocampal pyramidal cells (Avery and Johnston, 1996), a 
maximum value of 1.541 ms being reported for thalamic relay cells 
(McCormick and Huguenard, 1992). These values seem compati- 
ble with those given for the appropriate subtypes in Table 6. For 
skeletal muscle, a value of 19.82 ms (Standen and Stanfield, 1982; 
here assuming L-type) is also compatible with the relatively high 
value for the Cayl-l subtype. The experimental values reported for 
smooth muscle are around 4.0 ms (Muinuddin et al., 2004). Some 
available data for Tm for subtypes are given in Section 5. 

In Table A3-Modeling, there are reported 23 values employed 
for Tm, with a large range of values, many of which are much larger 
than the values for subtypes reported in Table 6. Many of the values 
employed for a given neuron type are not based on experiments for 
that neuron type but adapted from those for other, sometimes quite 
different, neurons. Often quite disparate values appear for what are 
possibly the same kind of cell. Discrepancies may arise because it is 
difficult to categorize a specific channel type in a neuron without a 
detailed voltage-clamp and pharmacological analysis. We consider 
cells in the following categories. 

Dopaminergic neurons. The value of about 0.134 ms was employed 
by (Li et al., 1996) which contrasts greatly with the value 19.5 ms 
used by both Amini et al. (1999) and Komendantov et al. (2004). 
The approach of Li et al. (1996) is critically discussed by Amini et 
al. (1999). 

Cardiac cells. Concerning the listed values, in the three modeling 
studies on ventricular myocytes, two (Luo and Rudy, 1994; Shan- 
non et al., 2004) have Tm,- values with a maximum around 2.25 ms 
whereas a third puts the maximum value much higher at about 19.3 
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ms (Fox et al., 2002) and at a considerably more hyperpolarized 
level. In their analysis of the firing properties of rabbit sino-atrial 
node pacemakers, Zhang et al. (2000) chose a value of about 3.8 ms 
for the maximum activation time constant. Thus in most of these 
studies the magnitudes of Tm are comparable with or near the values 
for Cavl.2 and Cayl.S channels (see Table 6) which are plentiful in 
the heart. 

Hippocampal cells. The range of activation time constants for hip- 
pocampal pyramidal cell models is from 0.83 ms for CAl (Poirazi 
et al, 2003) to 4.56 ms for CA3 (Migliore et al., 1995). For the 
latter cell-type, Jaffe et al. (1994) used a value with a maximum of 
1.5 ms. According to Table 6 and Figure 12 all of these values are 
within the ranges of the ones given for Cavl.2 and Cavl.3 channels, 
both of which are prevalent in the hippocampus. 

Thalamic relay cells. The experimental findings of McCormick and 
Huguenard (1992) for r^, with a maximum of 1.541 ms at l^ = —15 
mV, were implemented in a model by the same authors. The same 
values were later employed by Hutcheon et al. (1994). Rhodes and 
Llinas (2005) put the maximum value of the time constant some- 
what larger at 2.1 ms at -19 mV, but in the table of parameter 
values the current is referred as a generic Ica rather than specifi- 
cally as an L-type current. All of these values are compatible with 
the expected subtype values. 

Pituitary corticotrophs. The two modeling studies of these cells 
listed in the table have fairly large values of Tm- In the first study, 
LeBeau et al. (1997) set the value at 27 ms for all V, whereas 
Shorten and Wall (2000) had a maximum for r^ of value 11.32 oc- 
curring at V = —12 mV. There seems to be a fairly large difference 
between these values and the sub-type values, but there are many 
sources of variability as discussed elsewhere in this review. 

Cortical pyramidal cells. In the one parameter set available (Rhodes 
and Llinas, 2001), the maximal time constant for L-type activation 
was set at 2.1 ms, occurring at -19 mV. This value is well within the 
range of the available data for the appropriate subtypes as given in 
Table 6 and Figure 12. 

Spinal motoneurons. The L-type activation listed in the 3 listed 
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modeling studies of these neurons (Booth et al, 1997; Carhn et al. 
2000; Bui et al, 2006) is much slower than the values available for 
either of the subtypes Cavl.2 Cavl.3 which are candidates for being 
present in CNS neurons. The first of these studies set Tm = 40 
ms, whereas the two subsequent ones put r^ = 20 ms whereas the 
largest value in Figure 12 is about 4 to 5 ms. 

Smooth muscle. In modeling a rat mesenteric smooth muscle cell, 
Kapela et al. (2008) for Tm{V) used a Gaussian-like function with 
a maximum value of 3.65 ms at \^ = —40 mV. Such a value is close 
to the experimentally reported time constant in (Muinuddin et al., 
2004). 

Figure 5 shows two plots of Tm versus V for two of the few 
available experimental data sets (but see Section 5 for results for 
L-subtypes). These are both for rat sympathetic neurons (A, Bel- 
luzzi and Sacchi, 1991; B, Sala, 1991) and are, in the absence of 
reliable L-type data in physiological or near physiological solutions. 
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Figure 5: Activation time constants for high threshold Ca^^ cur- 
rents. Two sets of experimental results for sympathetic neurons. 
In both cases the corresponding steady state activation function is 
shown, using the same scale. A. From Belluzzi and Sacchi (1991) 
with 5mM Ca2+. B. From Sala (1991), with 4mM Ca2+. 
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Figure 6: Histogram of the maximum values of time constants r^, 
in ms, for activation for L-type Ca^"*" currents using data from Table 
A3. Includes both experimental and modeling studies less than 10 



ms. 



taken from HVA studies which may contain L-type and other high 
threshold components. These data sets are very different in mag- 
nitude and form from those used in many modeling studies, where 
sometimes the time constant for activation is taken to be a large 
value, constant and independent of voltage. They are also consid- 
erably less than the values reported for Cavl.2 and Cavl.3 in 2mM 
Ca^^ in Section 5. Molitor and Manis (1999) report experimental 
values of r^ with a maximum of about 1.2 ms at about V=-20 mV 
in guinea pig dorsal cochlear nucleus neurons. In Figure 5 are also 
drawn the steady state activation functions, also based on experi- 
ment, for comparison. 

The data of Table A3 on maximum values of activation time 
constants are collected into a histogram in Figure 6, but this in- 
cludes only values less than 10 ms. The mean for the experimental 
results on neurons is 1.34 ms. For all the data shown (n=17) the 
mean is 2.38 ms and the standard deviation is 1.51. 

For time constants Th for voltage dependent inactivation in neu- 
rons there are very few data. The average of the two available ex- 
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perimental maximum values is 65 ms and the average of 5 modeling 
maximum values is 140.9 ms. 

3.3 Magnitudes of IcaL or related quantities 

Data on the magnitudes of IcaL in various cells is collected in Table 
A5. Such data are presented in several forms, viz, (a) the actual 
current density, usually in pA/pF, which may be for the whole cell 
or a part thereof (b) the permeability (c) the conductance per unit 
area ql (d) the whole cell conductance Gl (d) the fraction of total 
calcium current that is IcaL- It can be seen that the proportion of 
IcaL across cell types is extremely variable from nearly all of the 
calcium current in some muscle and cardiac cells (Kapela et al., 
2008; Benitah et al., 2010) to just a few percent, for example in 
serotonergic (dissociated) cells of the dorsal raphe nucleus (Pening- 
ton et al., 1991). The roles of IcaL are presumably very different in 
these examples. From Table A5, the largest IcaL current density is 
19.3 pA/pF in supraoptic nucleus (Joux et al., 2001), but almost 
as high magnitudes are found in rat CAl pyramidal cells (Xiang et 
al., 2008) and ventricular cardiomyocytes in rats with renal failure 
(Donohoe et al, 2000). The highest reported conductance density 
is 7 mS per square cm in rat CAl pyramidal cells (Xiang et al., 
2008). Note that current densities can differ significantly in differ- 
ent preparations such as those with different ages, or with different 
times in culture. 



4 Examples of calculations of Icai 

In this Section we apply the basic model (1) to determine the L-type 
calcium current in three cases. In the first, only VDI is considered, 
whereas in the remaining two examples, VDI and GDI are both 
operative. 
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4.1 Steady state currents obtained by the lin- 
ear and constant field methods with voltage 
dependence only 

The activation and inactivation variables, which are dimensionless 
and take values in [0,1], may depend on membrane potential or ion 
concentrations, or both, but in this section only voltage dependence 
is considered. Here we wish to compare current-voltage relations 
obtained by the two methods outlined in Section 2.4. In order to 
compare the steady state currents IcaL,Lm and IcaL.CF calculated by 
the two methods, an average maximal conductance density g^ (from 
the values given in Table A5) of 1.3925 mS/cm^ was assumed for 
the determination of IcaL.Lin, although this value varies greatly from 
cell to cell. For the steady state activation variable rrioo used in the 
calculation of IcaL.Lin, two values of V^ i were employed, being a 
relatively low value, V^ i = —29.6 mV, and an intermediate value 
Vj^ 1 = —14.4 mV. The corresponding values of km were 6.97 mV 
(low) and 7.78 (intermediate). When inactivation was included, the 
values for V^^ i were \4 i = —43.5 mV (low) and Vf^ i = —36.2 mV 
(intermediate), with corresponding values for k^ of 9.55 mV (low) 
and 6.18 (intermediate). 

Four sets of values for pi and p2 were employed, giving the forms 
m,m'^,'mh and m'^h. The calculated results are shown in the four 
panels of Figure 7. In each case the current IcaL.Lin (linear method) 
was calculated first, the temperature being set at 25° C, the internal 
calcium concentration Cai fixed at 100 nM and the external concen- 
tration Cao fixed at 2 mM, giving Vca = 128.4 mV which was used 
as the reversal potential. The constant field current IcaL,CF was 
then calculated and the permeability adjusted to make the peak 
(inward, negative) current the same as that in the linear case. This 
gave 8 values of the permeability which are listed in Table 4 which 
are close to those reported for L-type calcium currents (Hutcheon 
et al., 1994; Luo and Rudy, 1994). 

The following can be seen from the results in Figure 7. Details 
of the values of V at which minima in IcaL (maximal current am- 
plitude) occur are given in Table 5. Firstly, without inactivation 
(top two panels), as expected, the constant field and linear calcula- 
tions are similar in magnitude until a value of V just greater than 
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that at which the minimum of IcaL (maximum amphtude) occurs. 
However, the minima in the constant field case occur at voltages 
about 11-13 mV more negative than for the linear case, both for 
open probabilities given by m or m^ and whether V^ i is ascribed 
the low or high value. For larger values of V the two calculations 
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Figure 7: Calculated asymptotic steady state (t = oo) currents 
for L-type Ca^^ currents with average activation and inactivation 
properties for a relatively low voltage activation (V^ i = —29.7 
mV), blue curves, and an intermediate value {V^ i = —14.4 mV), 
red curves. The solid curves are for the constant field current (for- 
mula (17)) and dashed lines are for the linear current case (formula 
(15)). The top two panels are for activation only (m and m^) and 
the bottom two panels are for activation with inactivation {mh and 
im?h). In each of the 8 sets of curves, the maximum amplitudes of 
the constant field and linear currents are scaled to be equal, for a 
conductance density of 1.3925 mS/cm^ in the linear case. Temper- 
ature 25°C. 
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Table 4: Estimated Pl from IcaL,CF in 


10 ^ cm/sec 


Kn,i 




m 


m' 


mh 


m^h 


Low (-29.6 mV) 
Intermediate (-14.39 


mV) 


3.01 

4.41 


3.40 
5.26 


2.08 
2.11 


2.41 
2.76 



give strikingly different results, as the constant field currents are 
almost symmetric about their minima and return to almost zero at 
values of V around 70 mV whereas the linear calculation has the 
current returning to zero very gently by comparison, and of course 
reaching zero at the reversal potential of Ca^+ (about 128 mV). 
Whether such differences are significant for any given cell depend 
on how long the voltage stays at values greater than those at the 
minima (maximum amplitude), which are given in Table 5. 

Changing from an m to an m^ dependence for IcaL shifts the 
voltage at which the greatest current amplitude occurs by about 5 
or 6 mV in the depolarizing direction. Furthermore, there is not a 
large difference in the magnitude of the current between the low and 
high V^ 1 cases. Secondly, with inactivation (bottom two panels) 



m,2 



Table 5: Membrane potentials in mV at maximal current amplitude 





K„,r 


=-29.7 mV 


Kn- 


=-14.4 mV 




Linear 


Const Field 


Linear 


Const Field 


m 


-9 


-20 


6.5 


-7 


m^ 


-4 


-15 


12 


-1 


mh 


-32 


-35 


-32 


-34 


m^h 


-25 


-27 


-18 


-21 



the current/voltage relations are practically symmetric for all com- 
binations of V^ 1 and the factors mh and m?h. Thus, the constant 

'"-'2 

field calculations are practically the same as the linear calculations 
for all values of V . Inactivation, either as mh or m^h shifts the 
voltage of maximal current amplitude by very large amounts in the 
hyperpolarizing direction. The smallest shift is -12 mV with m^h 
and the constant field calculation and the largest is -30 mV with 
m?'h and the linear calculation. The magnitudes of the currents are 
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reduced by an order of magnitude by the inactivation and the reduc- 
tions are much greater for the higher value of V^ i (red curves). The 
difference for the hnear and constant field results for the voltage of 
maximal current amplitude is only a few mV. 

4.2 Time dependent IcaL for high and low group 
neurons with VDI: m or m^? 

The question arises as to when the time-dependent results for an mh 
form of current differ significantly from those for an rn^h form. This 
was investigated using the parameters for steady state activation 
and inactivation for the two sets called high and low summarized 
in Table 3. 

The time constant of activation was set at the average maximal 
experimental value of 1.34 ms for neurons and the time constant 
for inactivation was set at 29.9 ms as found in Belluzzi and Sacchi 
(1991), this being one of the only two experimental values available. 
Current flows were computed for voltage steps from an initial value 
Vq = —80 mV to final values Vi = —30 mV and Vi = mV. The 
following Ho dgkin- Huxley (1952) formula may then be applied 

I{t;Vo,Vi) = gLiV^-Vre,,L)[m^-{m^-mo)e-'/--r 

where we have employed the abbreviations itlq = mao{VQ),mi = 
rriooiVi) for the activation, ho = /ioo(Vo))^i = ^oo(^i) for the inac- 
tivation, Tm and Th being the time constants. For these calculations 
'qi was set at unity as only relative magnitudes are of interest. 

The results are shown in Figure 8. It can be seen that in most 
of the examples shown, there are only minor differences in the cur- 
rents computed by the mh versus m^/i form. The case where the 
difference is very significant (plot 8A) is for the high neuron group 
with a voltage step from - 80 mV to -30 mV. However, only a de- 
tailed neuron model, which is beyond the scope of this article, could 
ascertain whether significant differences arise in the ongoing spiking 
activity when L-type currents are one of an array of current types. 
Nevertheless, it has been demonstrated that the choice of the form 
of the open probability contribution from voltage dependent terms, 
m^^h^'^, can have significant consequences for the current /voltage 

37 




-20 



l(t) 



50 100 150 200 





50 100 150 200 




50 100 150 200 50 100 150 200 

t in ms 



Figure 8: In all plots, current in arbitrary units versus time in ms. 
Blue curves are for p = 1 and red curves are for p = 2. A. High 
neuron group parameters (see Table 3) with a step from -80 mV 
to - 30 mV. B. High neuron group parameters with a voltage step 
from -80 mV to mV. C Low neuron group parameters with a step 
from -80 mV to -30 mV. D. Low neuron group parameters with a 
voltage step from -80 mV to mV. 

relations or for current versus time. In the ideal situation voltage- 
clamp data for both activation and inactivation would be used to 
find the best fitting theoretical form for a component current, as 
in the original Ho dgkin- Huxley (1952) analysis. Care is required 
in reporting parameters for activation and functions so that it is 
clear whether they are for m or m^. It is also important that the 
experimental data have been analyzed correctly. 
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4.3 Two time constants in IcaL inactivation 

For L-type calcium current inactivation, most reports indicate that 
there are two time constants, one fast and one slow (Luo and Rudy, 
1994; Romanin et al., 2000; Meuth et al., 2001; Budde et al., 2002; 
Lacinova and Hoffman, 2005; Faber et al., 2007; Findlay et al., 
2008). According to Lacinova and Hoffman (2005), such is the case 
for smooth muscle, ventricular myocytes, several neurons and the 
Cavl.2 channel. Furthermore, in cells, fast time constants range 
from 7 to 50 ms and slow ones from 65 to 400 ms whereas for the 
Cavl.2 channel, the ranges are 20 to 100 (fast) and 160 to 2000 
(slow). See Table A3 for several data although no distinction is 
made between fast and slow inactivation time constants in the cited 
works. In an earlier study of calcium currents (not specifically stated 
as L-type) in CAl pyramidal cells, Kay (1991) found, with 5mM 
Ca^"*", a fast time constant with a maximum value of about 240 ms 
at V=-10 mV, and a slow time constant with a maximum of about 
2200 ms at V=-ll mV. It was pointed out that typically (in cells 
without GDI) the amount of inactivation increases with increasing 
V, and that the rate of inactivation tends to be Gaussian-like as a 
function of V, with a maximum around the voltage V^ i at which 
the steady state inactivation is half-maximal. 

The relative contributions of GDI and VDI vary. Lacinova and 
Hoffman (2005) state that in general for a brief pulse, GDI is re- 
sponsible for about 80% of inactivation and VDI responsible for the 
remaining 20%. The faster time constant is usually taken to be 
associated with GDI. However, in some cases there are taken to be 
two time constants for VDI of L-type Ga currents, such as in ven- 
tricular myocytes (Faber et al., 2007; Findlay et al., 2008) where the 
dichotomy is attributed to the existence of two voltage-dependent 
inactivation states. Two time constants were also demonstrated by 
Morad and Soldatov (2005) for barium currents with no GDI for 
Gavl.2. 

To illustrate the occurrence of two time constants we show re- 
sults obtained for IcaL for the basic model of Equation (1) using 
the constant field expression for the current. These are designed to 
be comparable to the experimental results on thalamic relay (LGN) 
neurons in Budde et al. (2002), where a step to a test voltage of 10 
mV is made for 800 ms after holding at -60 mV. The parameters 
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Figure 9: Showing two time constants for calculated L-type Ca^"^ 
currents evoked by a test pulse of duration 800 ms at 10 mV after a 
holding potential of -60 mV. The blue curves are for GDI and VDI, 
whereas the red curves are for no GDI. The figure on right shows 
clearly the change in time constant from early to late (blue curve) 
compared with the single time constant for VDI only (red curve). 



used were as follows. For activation, VI, iwas set at -37.5 mV with 
slope factor 5 mV, and for volt age- dependent inactivation, V^ i was 
set at -59.5 mV with a slope factor of 10 mV. The time constant 
for activation was given by 



Tm = l + 2e 



-[{V-V 



i)/15p 
2 



(20) 



and that for VDI was set at Th = 300 ms for all V. In addition, a 
calcium pump term was added of the form 



Jpump\^"'i) ^c 



v^ai 

Gai + Kp 



(21) 
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where Kp = 0.0005 mM (Rhodes and Lhnas, 2005) and K^ = 
0.000002. Initial internal calcium concentration was set at 70 nM 
and extracellular at 2niM. The standard value of the multiplicative 
factor AplPI in (19) was chosen to make the maximum amplitude 
of IcaL about 0.2 nA as in Budde et al. Figure 2. For the GDI term 
we use the standard model (7) and (8) with n = 1, rj = 10 ms 
(constant) and Kf = 0.0002 mM. 

Results are shown for IcaL in Figure 9, the contribution from 
the pump being small and not included. During the pulse the in- 
ternal calcium concentration rose to a maximum of about 300 nM. 
In the left panel, IcaL is plotted against time in the cases of GDI 
+ VDI (blue curve) and for VDI only (red curve). With GDI the 
initial decline of current is very rapid and merges into a slower de- 
cay after about 100 ms; with VDI only, the decline in current is 
slow and relatively uniform as can be seen in the right hand panel 
where ln(— Ica,L) is plotted against time to reveal changes in time 
constant. The curve for VDI is a straight line, indicating a single 
time constant. The early (fast) time constant for GDI + VDI was 
approximately 49 ms and the slow time constant about 264 ms. 
Both of these values are in the ranges given by Lacinova and Hoff- 
man (2005) quoted above. They are also similar in magnitudes to 
the values reported in Budde et al. (2002). 

4.3.1 Effects of varying parameters 

The effects of varying several parameters on the inactivation of IcaL 
according to the above test pulse procedure are shown in Figure 
10. Here both GDI and VDI are present, but the graphical results 
are given only for the first 200 ms. The parameters used here are 
as given by Rhodes and Llinas (2005) as in Table Al.3, using the 
same /oo, but with time dependent /(Gai,t). The standard set of 
parameters is taken to be as for Figure 9, with the following varia- 
tions. In the top left panel, the current is calculated for four values 
of the GDI time constant rj, being 5, 10, 30 and 100 ms. Note that 
most modeling of L-type currents effectively takes tj = (see Table 
Al.3), as the steady state is assumed to be attained immediately. 
When Tf is 5 or 10 ms (blue and red curves), the current has almost 
the same time course as it does when the immediate steady state 
assumption is made. For r/ = 30 ms (green curve), which is the 
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value used by Standen and Stanfield (1982), the decay of current 
is much slower for the first 100 ms, very different from that for the 
immediate steady state assumption. For r/ = 100 ms (black curve), 
the decline does not reach values comparable with the remaining 
three cases for at least 200 ms. In the top right panel are shown 
IcaL for four values of the time constant r^ for voltage-dependent 
inactivation, no dependence on V being assumed. The value of Tf 
is set at 10 ms. Results are shown for time constants r^ (in ms) of 
100 (blue), 200 (red), 300 (green) and 400 (black), which values en- 
compass approximately the values employed (Table A3). As could 
be anticipated, increasing t^ from 100 to 400 has a negligible effect 
for the first 20 ms but quite noticeable effects on the late current. 
In the bottom left panel are shown results of varying the param- 
eter Kc which determines the rate of calcium pumping as in (21). 
The values employed are 1,2,3 and 10 times a base rate with blue, 
red, green and black curves respectively. When the pump rate is 
higher, Ca^^ is less available for GDI so the current has a larger 
magnitude. However, in the situation modeled here, it can be seen 
that a significant difference in IcaL occurs when the pump rate is 
increased tenfold, but that increases of 100% and 200% have only a 
small effect. Finally, in the bottom right panel are shown the time 
courses of IcaL for various channel densities, being 1 (blue), 4 (red), 
8 (green) and 12 (black) times a base rate. The current amplitude is 
assumed proportional to the channel density so the higher channel 
densities result in more effective GDI. These results indicate that 
these four quantities have a significant influence on the amounts of 
GDI and VDI and thus the magnitude and time course of IcaL- 

4.3.2 CICR 

If there is a substantial amount of GIGR, as occurs in cardiac and 
other muscle cells, then there may be two phases of GDI, consist- 
ing of an early rapid component due to release of Ga^+ from the 
sarcoplasmic reticulum and then a slower phase in response to cal- 
cium ion entry through L-type channels (or other calcium channels) 
(Budde et al., 2002; Bodi et al, 2005; Gueni et al., 2009; Empson 
et al., 2010). This raises the possibility of three time constants 
in the decay of IcaL, two for the two phases of GDI and one for 
VDI, or even four time constants if as mentioned above there are 
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Figure 10: The effects of changing 4 parameters on the time course 
of the L-type calcium current when there are both GDI and VDI. 
Top left: the time constant Tf for GDI is varied. Values in ms are 
5 (blue), 10 (red), 30 (green) and 100 (black). Top right the time 
constant Th of VDI is varied. Values in ms are 100 (blue), 200 (red), 
300 (green) and 400 (black). Bottom left: The effects of various 
calcium pump strengths K^ on IcaL - see Equation (21). Values of 
Kc are 1 (blue), 2 (red), 3 (green) and 10 (black) times a base rate of 
0.000001. Bottom right: How the channel density p^, through the 
magnitude of the current, affects IcaL- The multiplicative factor 
ApiPl has values 0.25 (blue), 1 (red), 2 (green) and 3 (black) times 
the base rate, where the base rate is the standard value described 
for the previous figure. 

two for VDI. The effect of Ga^+ released from sarcoplasmic retic- 
ulum is graphically illustrated in Hinch et al. (2004). The basic 
model (1) can incorporate the influence of GIGR on GDI by tak- 
ing account of the time-course of the Ga^"*" concentration, as stated 
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in Section 2.1. This quantity will usually be obtained from either 
ordinary (for Cai(t)) or partial differential equations (for Cai(x, t), 
where x is from 1 to 3 space variables) which incorporate calcium 
ion concentration changes due to pumping and buffering as well as 
fluxes through various channels. See Luo and Rudy (1994) for such 
a deterministic example. A detailed stochastic model for CICR in 
cardiac myocytes has been analyzed by Williams et al. (2007). 

4.4 Calculations for double-pulse protocols 

According to Budde et al. (2002), GDI tends to be rapid and results 
in a "U-shaped" inactivation curve when investigated using a double 
pulse protocol (Tillotson, 1979; Kay, 1991; Meuth et al., 2001). 
We show corresponding results for the standard model for L-type 
calcium currents as described in section 2. The voltage is held at a 
"conditioning" value, Vcond followed by a pause and then a "test" 
pulse at Vtest is delivered. Of interest is the relationship between 
the first and second responses. As in Figure 1 of Budde et al. (2002), 
mathematically we have 



Vit) 



Vi, < t < ti, 

Vcond, ti <t < t2, 

Vi t2<t< ts, (22) 

Vtest h <t < ^4, 

Vi t > ti. 



We use the constant field method and assume, as the results of 
Section 4.2 indicate a relative insensitivity to the power of m, 

IcaL = F(\/, Cai, Cao)m2(\/, t)h{V, t)f{Cs.u t) (23) 

where F contains physical constants and the constant field factor as 
given by (17). For the GDI term we use (7) and (8) with n = 1, 
Tf = 30 mS (constant) and Kf = 0.0005 mM. 

For the volt age- dependent activation and inactivation parame- 
ters we use the two sets for thalamic relay cells as given by Pospis- 
chil et al. (2008) (see Table Al.l) and Rhodes and Llinas (2005) 
(see Table Al.3). The time constant for activation was chosen as 
in Equation (20), and the inactivation time constant was set at 
Th = 500 ms for all V. 
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Figure 11: Calculated L-type calcium currents in the two-pulse pro- 
tocol described in the text using the standard model given by (1) 
with VDI and GDI. Time is in ms, voltages in mV and currents 
in nA. In the left column are results with data for activation and 
VDI of Pospischil et al. (2008), and in the right column of Rhodes 
and Llinas (2005). The voltage is held at Vcond from t = 20 to 
t = 220 and at Vtest from t = 225 to t = 445. In the bottom 
graphs the amplitude (minimum) of IcaL is plotted against the con- 
ditioning pulse voltage for the first pulse (blue curves) and for the 
second pulse (red curves). Vtest is at the minmum of the blue 
curves which occur at -20 mV for the left set of results and -5 mV 
for the right set. In the top graphs, responses are shown for three 
particular values of Vcond- The red curves are for Vcond at the 
value of Vtest, which gives a large reduction in IcaL on the second 
pulse due to GDI. The blue curves are for Vcond = 70 where the 
current on the first pulse is very small. The green curves are for 
Vcond = —35 (left) and Vcond = —25 (right). The results may be 
compared with those of Figure 1 in Budde et al. (2002). 
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The conditioning pulse is on from 20 to 220 ms and the test 
pulse is on from 225 to 445 ms. The voltage at the test pulse is 
taken as that at which the response to the first pulse has the largest 
(negative) amplitude. Results of these calculations are shown in 
Figure 11. The curves in the bottom part of the figure show the 
magnitudes of the currents in response to the first (blue) and second 
(red) pulses. The reduction in the amplitude of in the second pulse 
is most prominent when the response to the first pulse is the great- 
est, because Cai has increased the most, resulting in the greatest 
amount of GDI. In the top parts of the Figure are shown responses 
to two-pulse clamps for three particular values of Vcond, as ex- 
plained in the Figure caption. Although such results indicate that 
the standard model for IcaL with the usual voltage-dependent acti- 
vation, coupled with VDI and GDI may be adequate for describing 
certain experimental paradigms, it is uncertain whether this is the 
case in applications to single cell models, mainly due to the lack of 
concrete data. The relative contributions of GDI and VDI depend 
on the particular cell type, mainly through the various combina- 
tions of subtypes of L-channels as well as modulating factors such 
as sub-unit configurations and various ligands. 

4.5 Difficulties in the accurate modeling of IcaL 

Mathematical modelers of nerve cell physiology are invariably faced 
with the dilemma of insufficient accurate data concerning the many 
ion channels influencing a neuron or muscle cell's activity. It is 
hard to think of a single neuron for which such data is available 
for all known channel types. As an example, L-type current data, 
for a current denoted by II, were not available for thalamic relay 
neurons when they were first comprehensively modeled (McGormick 
and Huguenard, 1992). Instead, the parameters were obtained from 
the calcium current results of Kay and Wong (1987) for GAl pyra- 
midal cells, where in fact an L-type current had not been pharmaco- 
logically identified. Kay and Wong hypothesized that the channels 
they had analyzed were from a homogeneous population. The same 
or similar data were employed by Traub et al. (1991) for Ica in 
GAl pyramids, Hutcheon et al. (1994) for I^ in thalamic relay 
cells, Traub et al. (1994) for for Ica in GAS pyramids, Wallen- 
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stein (1994) for J^ in nucleus reticularis tlialami neurons and Wang 
(1998) for Ica in cortical pyramidal cells. In none of these models 
was calcium-dependent inactivation included. 

It is hoped that in the future, experiments can be performed to 
remedy the paucity of accurate data required for L-type (and other) 
current modeling. There are other limitations apart from the use of 
data from one cell type for another. The fact that the majority of 
experiments are performed on dissociated cells means that existing 
data must be biased towards those for cell somas, so that guesses 
have to be made for dendrites. To make choice of parameters even 
more difficult is the use of various different concentrations of either 
calcium or barium ions which may give results for such quantities as 
activation potentials which are very different from those for phys- 
iological concentrations. Temperature is another key factor which 
may alter the dynamic properties considerably. Furthermore, re- 
sults for heterologous ion channels may not be the same as for the 
corresponding channels in native cells. 

Thus, given the many uncertainties in the parameters describing 
the kinetics of the various ion channels, it is hard to know whether 
model-based predictions are truly descriptive of cell electrophysiol- 
ogy even when they agree with experimental observations. The hope 
is then that the modeling results are robust enough with respect to 
relatively small changes in the parameters. To construct a model 
of a nerve cell is a very large undertaking, and fortunately many 
researchers have had the courage to proceed with such endeavours 
in the light of uncertainties concerning the underlying channel data. 

5 L-channel subtype properties 

L-type currents are not only important for neuronal and cardiac 
cell dynamics. Although L-type channels play a limited role in 
the process of synaptic transmission, they are crucial for activity- 
dependent gene expression and for regulating plasticity at certain 
synapses (Hardingham et al. 2001; Dolmetsch et al., 2001; Helton 
et al. 2005). Thus they have been found to play an essential role in 
long-term alterations in synaptic efficacy underlying learning and 
memory in the hippocampus (Kapur et al., 1998; Graef et al., 1999; 
Leitch et al., 2009) where both Cavl.2 and Cavl.3 channel subtypes 
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are predominantly located in postsynaptic dendritic processes and 
somata. Furthermore, current through Cavl.3 channels has been 
demonstrated at the output synapses of mice All amacrine cells 
(Habermann et al., 2003). 

5.1 Distribution 

The Cavl.l subtype is mainly found in skeletal muscle and Cavl.4 
is found in the retina, but also in human T lymphocytes (Kotturi 
and Jefferies, 2005). The remaining subtypes Cavl.2 and Cavl.3 
constitute the main L-type channels in neurons, endocrine cells and 
heart cells. The Cavl.2 subtype is widely expressed in heart, brain, 
smooth muscle and endocrine systems (Ertel et al., 2000). In the 
brain it is found in cortex, hippocampus, thalamus, hypothalamus, 
caudate putamen and amygdala (Splawski et al., 2004). In neu- 
roendocrine cells, Cavl.2 and Cavl.3 channels are involved in ac- 
tion potential generation, bursting activity and hormone secretion 
(Lipscombe et al., 2004; Marcantoni et al., 2007). According to 
Vignali et al. (2006), 60% of calcium current in both mouse pan- 
creatic A- and B-cells is L-type, with mainly Cavl.2 in B-cells and 
both Cavl.2 and Cavl.3 in A-cells. The Cavl.3 subtype is found 
in sinoatrial node, cochlear hair cells, and dendritic neuronal pro- 
cesses (Dolphin, 2009). Schlick et al. (2010) reported distributions 
of subtypes in many brain regions, including during development, 
and found that patterns tended to be intrinsic rather than depen- 
dent on neural activity. Using barium as charge carrier, Navedo 
et al. (2007) investigated isoform contributions in mouse arterial 
smooth muscle. Both Cavl.2 and Cavl.3 were involved in calcium 
sparklet formation but Cavl.2 was primarily responsible (see Table 
A6). 

Immunofluorescence studies have determined that both Cavl.2 
and Cavl.3 have mainly proximal locations on neurons, particu- 
larly on the dendrites, and they may make a substantial contribu- 
tion to the inward current and the action potential (Holmgaard et 
al., 2008). According to Martinez-Gomez and Lopez-Garcia (2007), 
in mouse and rat spinal neurons Cavl.2 channels tend to localize 
in soma and proximal dendrites whereas Cavl.3 channels are also 
found in distal dendrites. However, Zhang et al. (2006) found that 
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in the cat, Cay 1.3 channels were dense in ventral horn motoneurons, 
occurring mainly on somata and proximal dendrites and being re- 
sponsible for plateau potentials in these cells. Forti and Pietrobon 
(1993) found that functionally different L-type Ca^^ channels co- 
exist in rat cerebellar granule cells. Although Cavl.2 and Cay 1.3 
are often found in the same neuronal processes, particularly den- 
drites, their subcellular distributions are sometimes distinct (Hell 
et al. 1993). An interesting graphical representation of the distri- 
bution of various calcium channels over the surface of hippocampal 
pyramidal cells was given in Magee and Johnston (1995). 

Differential effects of corticosterone have been found for Cavl.2 
and Cavl.3 in the CAl region of hippocampus and the basolateral 
amygdala (Liebmann et al., 2008). In the former, corticosterone 
increases the amplitude of the slow afterhyperpolarization whereas 
in the latter no such effect is observed. The different responses are 
thought to reflect the expression of Cavl.3 in hippocampus and the 
absence of this subtype in the basolateral amygdala. 

5.2 Role in pacemaking 

In the heart, Cavl.3 is found in atrial tissue where it contributes 
to pacemaking but not in ventricular muscle that expresses Cavl.2 
(Lipscombe et al. 2004). According to Torrente et al. (2011), SA 
node pacemaker activity controls heart rate and local Ca^"*" release 
is tightly controlled by Cavl.3 channels. In Cavl.3 KO mice, the 
frequency of Ca^+ transients in pacemaker cells was reduced by 45% 
relative to WT. In another study with KO mice, Zhang et al. (2011) 
found that Cavl.3 channels are highly expressed in atrioventricular 
node cells and played a critical role in their firing. Mesirca et al. 
(2010) pointed out that atrioventricular node cells may contribute 
to pacemaker activity in the case of SA node failure. 

The properties of Cavl.3 channels which make them more suit- 
able for pacemaking activity than Cavl.2 channels were discussed 
in Vandael et al. (2010) with particular attention to dopaminer- 
gic neurons, cardiac pacemaker cells and adrenal chromaffin cells. 
Cavl.3 activates with steep voltage dependence at voltages con- 
siderably more negative than Cavl.2 and has faster activation but 
slower and less complete voltage-dependent inactivation. According 
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to Striessnig (2007), such activation at low voltages enables Ca^LS 
Ca^"^ currents to support pacemaking activity in sinoatrial node and 
some neurons (for example Putzier et al., 2009a). For adrenal chro- 
maffin cells, Marcantoni et al. (2010) found that although Cavl.2 
and Cavl.3 occurred with about equal frequency in WT cells, most 
of which ffied at about 1.5 Hz, Ca^LS-deficient cells had signifi- 
cantly less ffiing. The essential role of Ca^LS and their coupling 
to BK channels in the pacemaker activity of these cells was estab- 
lished. Alternative splicing of the C terminus of the ai subunit of 
Cavl.3 has been shown to give two forms, of which one (short, 42 A) 
is expected to be more supportive of pacemaker activity (Singh et 
al., 2008) with 73.6% inactivation of Ica in 30 ms. See also Table 
A6. 

5.3 Biophysical properties: activation and VDI 

In Table A6 are listed some of the available data on voltage-dependent 
activation and inactivation for examples of the four main subtypes. 
As pointed out in the Introduction, there are no definitive results 
for each subtype because of the modulatory effects of other subunits 
and varying experimental conditions such as the use of different 
hosts and different ionic concentrations. That is, one can expect 
considerable variability in channel properties for a given subtype 
from different preparations. As seen in the examples shown in Fig- 
ure 3, Cavl.3 channels can have quite low activation thresholds 
compared to the Cavl.2 subtype and which are, in the data given 
here, similar to that of T-type. The latter have strong VDI whereas 
Cavl.3 currents are reported to have have weak VDI and strong 
CDI (Lipscombe, 2002). Very few modeling studies have distin- 
guished subtypes, an exception being Shapiro and Lee (2007). It is 
also apparent from the steady state inactivation function for Cavl.2 
shown in Figure 3, based on Hu and Marban (1998), that virtually 
no Cavl.2 channels would be available to open around resting mem- 
brane potential (about -60 mV). 

In a pioneering study using tsA201 cells, Koschak et al. (2001) 
determined electrophysiological properties for Cavl.2 and Cavl.3 
(see Table A6) with current densities of 24.5 and 24.3 pA/pF, re- 
spectively, in 15 mM Ba^"*". Generally in the cases reported Cavl.3 
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has faster kinetics than Cavl.2 (Hehon et al., 2005) which is seen in 
the data of Table 6 and Figure 12. Inspection of Table A6 reveals 
surprisingly that V^ i is lower in the example for Ca^LS than for 
Cavl.2, but that V;^ i is higher for Cavl.3 than for Cavl.2 . Unfor- 
tunately the incompleteness of the data for all subtypes prevents a 
quantitative comparison of their relative current amplitudes, either 
singly or in combinations. Such a program will be useful for future 
work when more data are available. 

Table 6: Examples of data on activation time constants for L- 
current subtypes 

Reference Cayl.l Cavl.2 Ca^LS Cay 1.4 

Koschak (2001) 2.3 ms 1.2 ms 

15mM Ba++ 15mM Ba++ 

Hildcbrand (2004) Control 1.34 UlS 

2mM Ba^^ 10 /iM allcthrin 1.23 UlS 

Catteraii (2005) > 50 ms 1 ms < 1 ms <1 ms 

10 mV 10 mV 10 mV 

Helton (2005) 2.43e-0-°2^ 1.26e-02^ 

Kg [-20, 20] 1/ €[-55,10] 

Shapiro (2007) Model, MN 30 mS 24 mS 

Luin (2008) young 53.6* 

old 62.5* 



There are only limited data on activation time constants for L- 
channel subtypes, some being given in Table 6. Those for Cavl.2 
and Cavl.3 are graphed in Figure 12. Data from Koschak et al. 
(2001) are for tsA201 cells. The data of Helton et al. (2005) were 
given over the membrane potential ranges —55 < Vm < 10 for 
Cavl.3 and — 20 < Vm < 20 for Cavl.2. The following expression 
was employed by Putzier et al. (2009a) for the activation time 
constant of Cavl.3: 

_ / 0.020876(V- + 39.726) ^ ^g^^^^-(v^,5.33s)/224.2iy^ ^24) 

^m-\ ]^ _ g-(y+39.726)/4.711 +U.iy^^^e I . [Z^) 

This function is plotted in Figure 12 and can be seen to be 
approximately the same as the function 1.26e~°^^ given by Helton 
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Figure 12: Solid lines: time constants for activation, in ms, versus 
membrane potential in mV, found for examples of L-channel sub- 
types Cavl.2 and Cavl.3. Red and blue curves from Helton et al. 
(2005) (see also Table 6), black curve from Putzier et al. (2009a) 
(see also Equation (24)). Dashed lines: corresponding activation 
functions rrioo, scaled upwards. 
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et al. (2005) over the range [—55, 10]. Also shown for comparison 
in Figure 12 are the steady state activations rriooiV) corresponding 
to the three graphs of r^ versus V. 

5.4 Pathologies linked to Cavl.2 and CaylS chan- 
nels 

Several pathologies have been linked to Cavl.2 and Ca^LS channels. 
Examples include autism, bipolar disorder, and Parkinson's disease 
(Andrade et al., 2009). On the other hand, Cavl.2 and Cavl.3 
were shown to not be directly related to absence epileptogenesis 
(Weiergraber et al., 2010). Mutations of Cavl.2 have been found to 
result in Timothy syndrome (Splawski et al., 2004), which is a multi- 
organ disease manifesting in cardiac arrythmias, metabolic and im- 
mune system deficiencies and autism. It has been posited (Barrett 
and Tsien, 2008; Yarotskyy et al., 2009) that strongly decreased 
voltage-dependent inactivation, but not CDI, associated with the 
mutant channel subtype is a possible underlying cause. Alterations 
in L-channel density and activity have been found in aged rat hip- 
pocampal neurons with the possibility of a role in cognitive decline 
and the occurrence of Alzheimer's disease (Thibault and Landfield, 
1996; Thibault et al., 2001, 2007; Norris et al., 2010). The expres- 
sion of amyloid precursor protein was shown to increase L-type Ca^+ 
currents, thus leading to the inhibition of calcium oscillations (San- 
tos et al., 2009). Further, performance of rats in a maze task was 
found to be strikingly negatively correlated with L-channel density 
(Thibault and Landfield, 1996). 

Altered properties of L-type calcium channels have also been 
implicated in muscular dystrophy (Collet et al., 2003) and epilepsy 
(Fuller-Bicer et al., 2009). N'Gouemo et al. (2009) found upregu- 
lation of various high threshold calcium currents including IcaL in 
rat models of epilepsy. 

The role of Cavl.3 in pacemaking in sino-atrial node and in- 
ner hair cells of the cochlea has been demonstrated in KO mice 
which exhibit significant sinus bradycardia and hearing loss (see 
Lipscombe et al., 2004). Striessnig et al. (2004) reviewed the role 
of L-type calcium channels in human disease caused by genetic chan- 
nel defects (channelopathies) in ai subunits. Recently it has been 
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demonstrated in humans that a mutation in a gene which encodes 
the tti-subunit of Ca^LS channels is associated with deafness and 
sino-atrial node dysfunction (Baig et al., 2010). The possibihty 
and usefulness of subtype-specific therapeutic drugs has been dis- 
cussed in Striessnig et al. (2006) and more recently in Zuccotti et al. 
(2011). A comprehensive review of L-subtype properties (and other 
calcium currents) and associated channelopathies was compiled by 
Piedras-Renteria et al. (2007). 

6 Concluding remarks 

L-type Ca^+ currents have many fundamental roles, including pro- 
viding basic depolarization in action potential generation, amplifi- 
cation of synaptic inputs (established in spinal motoneurons) and 
involvement in signaling pathways for activating transcription fac- 
tors such as CREB and hence the expression of genes that are essen- 
tial for synaptic plasticity and other important cellular processes. 
Cavl.4 L-type channels are also involved in Ca^^-signaling in im- 
mune cells 

Quantitative descriptions of L-type Ca^+currents thus play an 
important role in understanding the dynamics of many processes 
in neurons and other cells, particularly cardiac myocytes. We have 
here reviewed the basic deterministic methods for a mathematical 
description of these currents. Whereas activation is usually taken 
to be voltage dependent, inactivation has varying degrees of GDI 
and VDI and each cell type requires a separate analysis of their 
relative contributions, which requires experimental investigation in 
addition to fundamental current-voltage relations. Incorporating L- 
type Ca^+currents in a neuronal model is often hampered by lack 
of knowledge as to the location of the corresponding ion channels, 
which is important to ascertain their contributions to cell dynamics. 
An additional complication is that some cells have more than 1 
subtype of L-type channels with separate properties and locations. 

The use of Markov chain theory, not described here, for transi- 
tions between various channel states has proven useful for cardiac 
cells but as far as can be discerned has not yet been applied to 
L-type channels in neurons. The deterministic, or averaging ap- 
proach, outlined in this review has the advantage of simplicity and 
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is thus likely to continue to be employed in modeling L-type (and 
other) calcium currents in multi-component nerve cell models. It is 
nevertheless not capable of capturing the fine details of the gating 
schemes for the activation and inactivation processes (Peterson et 
al., 1999; Erickson et al., 2003; Tanskanen et al., 2005; Tadross and 
Yue, 2010; Williams et al., 2010) which are made complicated by 
the roles of subunits and the varying degrees of voltage-dependent 
and calcium-dependent dynamics. See also Section 4.4 for a discus- 
sion of some of the general difficulties encountered in constructing 
accurate mathematical models for neurons. 
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8 Appendix: Tables of data on L-type 
currents 

The following tables are representative of experimental and mod- 
eling studies involving L-type calcium currents or calcium currents 
with properties similar to (or claimed to be similar) to L-type. In 
the years following the introduction of the term L-type, the term 
was often taken to imply a high threshold current, although later 
some L-type currents were found to be activated at quite low mem- 
brane potentials. The pharmacological distinction between the var- 
ious types is desirable but has not always been possible. Except 
in the case of experiments on subtypes which are clearly identified, 
usually in transfected cells, original sources may be consulted to 
determine the likelihood that the ion channels under consideration 
belong to the L-type category. Table Al.l contains activation and 
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inactivation data and modeling studies as opposed to experimental 
determinations have been marked M. Table A3 on time constants 
has been divided into an experimental part and a modeling part. 
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Table Al.l 

Parameters for L-type Ca^"*" activation and, if included, VDI 

Non cardiac cells *Denotes estimated M denotes model 



Source, cell type 



Form for current V^ i 
if given 



"'2 



Fox (1987) 

Chick DRG, lOmM Ca^"*" 

Kay (1987) 

CAl pyramidal, 2mM Ca^" 

Fisher (1990) 

CAS pyramidal, lOOmM Ba 

Belluzzi (1991) 

SYMP N, 5mM Ca^+ 

Sala (1991) 

SYMP N, 4mM Ca^+ 

Thompson (1991) 

CAl pyramidal, 3-lOmM Ba 

Traub (1991), M 

CAS pyramidal, Vca=^*0 

McCormick (1992), M 

TR cell, 2mM Ca^+ 

Schild (1993), M^ 

Rat med, 2.4mM Ca^"*" 

Hutcheon (1994), M 

TR med dors, 2mM Ca^+ 

Jaffe (1994), M 

CAS pyramidal 

Traub (1994), M 

CAS pyramidal, 2mM Ca^+ 

Wallenstein (1994), M 

NRT N, 2mM Ca^"*" 

Johnston (1995) 

General 

Kuryshcv (1995) 

Rat corticotropcs, lOmM Ba' 

Migliore (1995), M 

CAS pyramidal, 2mM Ca^ 



2+ 



2+ 



m 



-19 



4 -40 

7.73 



17 4.7 

m^h, const, field -8.1 9.85 -19.91 4.5 



m 



m^h{V - Vca) 



-4 



-19 



m^h{V - Vca) 



9.7 



7.73 



-78 9.9 



m^, const, field -14.1* 8.7* 



m^hiV - Vca) 

m^, const, field -14 

m?, const, field 4 

mHV - Vca) 

m?, const, field -14.1 8.7 

rn^h, const, field -15 



2.8 9.85 -14.61 4.5 

9.1 

4.7 
-19 7.73 



-12.3 7.8 
-11.3 5.7 
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Table A 1.1 continued 



Source, cell type 


Form for current 
if given 


y,,.k 


KfYi 


^M 


kh 


Magee (1995) 




9 


6 






CAl pyramidal, 20mM Ba^+ 












Avery (1996) 




-30 


6 






CAS pyramidal, 2 mM Ca^+ 












Li (1996), M 


mHV - Vca) 


-20.0 


5.3 






DA SN, Vrev,L=^'20 












LeBeau (1997), M 


rn^, const, field 


-12 


12 






PC 












Booth (1997), M 


miV - Vca) 


-40 


7 






Spinal, Vrev,L=80 












Wang (1998), M 


m{V - Vca) 


-20 


9 






CORT pyramid, Vc'a=120 












Molitor (1999) 




-10 


8 






DCN N, lOmM Ba^+ 












Athanasiades (2000), M 


m^h{V - Vca) 


-22.8 


9.85 


-34.61 


4.5 


MR, 2.4 mM Ca2+ 












Guyot(2000) 




-15 


6 


-33.5 


4.5 


2 1 
Bovine adrenal, 20mM Ba ^ 












Shorten (2000), M 


m^, const, field 


-18 


12 






PC, 20 mM Ca^+ 












Carlin (2000), M (HVA) 


miV - Vca) 


-10 


6 






Carlin (2000), M (LVA) 


m{V - Vca) 


-30 


6 






Spinal MN, Vrev,L=^0 












Joux (2001) 




-27.3 








SON, 5mM Ca^+ (Double Bolt) 




-12.4 








Zlmravleva (2001) 




-43.4 


5.0 


-78.2 


11.5 


LD thalamic N, lOmM Ba^+ 












Collet (2003) 


m{V - Vca) 


-1.1 


6.7 


control 




Mouse, SM, 2.5 mM Ca^+ 




-6.4 


6.2 


mdx 




Schnee (2003) 




-35 


4.7 


high frequ 




Turtle AHC, 2.8 mM Ca^+ 




-43 


4.2 


low frequ 




Durante (2004) 




-31.1 


5.5* 






Rat DA SNe, 2mM Ca^+ 












Muinuddin (2004) 


distal 


-7.8 


16.8 


-33.4 


11.6 


Cat SMM, 20mM Ba^+ 


proximal 


-6.1 


17.3 


-36.5 


13.2 


Jackson (2004) 




-36 


5.1 






Dorsomedial SCN, 1.2 mM Ca^+ 
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Table A 1.1 continued 



Source, cell type 



Form for current V^ 
if given 



"'2 



Liu (2005) 

Rat IC, 2 mM Ca^+ 

Bui (2006), M 

Spinal MN, Vrev,L=GO 

Marcantoni (2007) 

Adrenal CH cells, 2 niM Ca^+ 

Kager (2007), M 

CAS pyramidal, 1.5 niM Ca^+ 

Kapcla (2008), M 

Mesenteric SMM, 2.5 mM Ca^"*" 

Luin (2008) 

Cultured human SM, 10 mM Ca^+ 

Pospischil (2008), M 

TR, CORT, Fre«,L = 120 

Xiang (2008) 

CAl pyramidal, 5 mM Ca^"^ 

Marcantoni (2009) 

Adrenal CH cells, 2 niM Ca^+ 

Marcantoni (2010) 

Adrenal CH cells, 2 mM Ca^+ 

Adrenal CH cells, 2 mM Ca^"*" 



-16.6 



m 


-33 




-30* 


m^h{V - Vca) 


-11.3 


mh const, field 


-40t 




Printed as 


young 


8.73 


old 


14.47 


m'^h{V - Vca) 


-33.4 



-15 



2.45 -16.8 5.95 

6 

9* 



5.7 



8.3 



-42 



9.1 



7.72 -5.77 7.30 

7.74 2.38 10.1 

4.45 -58.5 20.28 

16.7 -30.5 11.2 

11* 



-23 (WT) 6.9 
-18.1(K0) 7.9 
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Table Al.2 
Parameters for L-type Ca^"*" activation and, if included, VDI 

Cardiac cells 



Source, cell type 


Form for current 
if given 


^m,l 


•^m 


^M 


kh 


Lindblad (1996), M 


{mih + 777-2) X 


-0.95(1) 


6.6(1) 


-28 


4.52 


Rabbit AM, 2.5 mM Ca^+ 


x{V-Vca) 


-33 (2) 


12(2) 






Handrock (1998) 




-3.5 


4.9 


-60 


7.5 


Failing human VM, 1 mM Ca 












Van Wagoner (1999) 




-14* 


7.75* 


-31.7 


4.6 


Human AM, 1 mM Ca^"*" 












Donohoc (2000) 








-25.9 




Rat VM, 1 mM Ca^+ 








-25.1 RF 




Zhang (2000), M 


[mih + 7712) X 


-23.1 (1) 


6 


-45 


5 


Rabbit SA PM, 2 mM Ca^"*" 


AV-Vca,L) 


-14.1 (2) 


6 






Arikawa (2002) 


(Control) 


-10.7 


6.2 


-28.2 


6.2 


Rat VM 


(Diabetic) 


-11.5 


6.0 


-27.2 


6.4 


Mangoni (2003) 




-25 




-45 


10 


Mouse SA PM, 4 mM Ca^+ 












Mangoni (2003) 




-3 




-36 


4.6 


Mouse SA PM 1.3KO 












Benitah (2010) 




-15 


8.5* 


-35 


5.95 * 


VM, 1.5 mM Ca^+ 
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Table A1.3 
Parameters for L-type Ca^+ current activation, VDI and/or GDI. 

Neurons and secretory cells 



Source, cell type 



Form for current 1/ i 



"'TTt Vfi 1 and/or kfi 

equation for J (Gaj, t j 



Amini (1999), M 

MID DA, 2.4 mM Ca^+ 

Rhodes (2001), M 

CORT pyramidal 

Poirazi (2003), M 

CAl pyramidal soma 

Komendantov (2004) 

MID DA, 2.4 mM Ca^+ 

Rhodes (2005), M 

TR, Vrev,L=60 

Komendantov (2007), M 

Mag Endocrine, 2.4mM Ca^"*" 



mfooiV - Vca) -50 



m- 



''hfoo{V-Vca) -19.3 t 



fnfoo const, field -30.1 
mf^iV - Vca) -50 



m 



'hfoo{V-Vca) -19t 



"i(/l,oo + /2,oo)x -27 
xiV-Vca) 



3 

7.5 

4.8 

20 

8 

4.5 



/c 



0.00045 



/< 



°° (0.00045+Cai) 

42 

0.0005 



oo 



f'. 



0.0005+Cai 

.001 
.001+Cai 

_ 0.00045 
°° (0.00045+Cai) 

42 

0.0001 
oo — 0.0001+Cai 

_ 0.2x0.0001"^ 

" O.OOOl^+Ca;" 

0.8x0.002 



U 



/c 



/c 
/l,oo 



2.00 0.002+Cai 
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Table Al.4 
Parameters for L-type Ca^+ current activation, VDI and/or GDI 

Cardiac and skeletal muscle cells 



Source, cell type 



Form for 
current 



VI 



km f^fi i and/or 



T7 

m J 

const, field 



equation for J (Cai, t J 



kh 



Standen (1982), M 
Insect SM, Ica 



Luo (1994), M 
VM, 1.8 niM Ca^"*" 

Fox (2002), M 
VM, 2 mM Ca^+ 



Shiferaw (2003), M 
VM, 1.8 niM Ca^+ 
Shannon (2004), M 
VM, 1.8 niM Ca^+ 

Stewart (2009), M 
PF, 2 niM Ca^+ 



-21.1 4.69 d/M = (/oo - /)/r/ 
/c - ' 



mhfoo 
const, field 

mf 

const, field 



mhfoo 
const, field 
mhf 
const, field 

mhih2f 
const, field 

See source for details 



-10 



-10 



6.24 



6.24 



5 6.24 

-14.5 6 

-8 7.5 



°o l+(Cac/X/) 

Cac= shell cone Ca^^ 



1 



^fO+Wf) 



Uf = 0.1ms ^ 
Kf = 1.0/iM 
« -35* 



l+(Cai/0.0006J 

-12.5 

df/dt = (/oo - f)/Tf 
J°° ~ T+(Ca~/Kjy 

Tf = 30, Kf = 0.18^M 
-35 

f _ .0005 
■1°° .0005+Cai 

-34.9* 

df/dt = 11.9(1 - /)- 

-1.7Cai/ 

-20 

/oo = 0.4 + , °-^ .. 

l+(Cai/0.05j 
80 



y = 2 + 



l+(cai/0.05)' 



8.6 



8.6 
3.7 
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Table A2 
Forward and backward rate functions, a and /3 for tti, and, if given, h 

All cell types 
Source, cell type am,h{V), (3m,hi^) '^rn n 



Standen (1982) am 

Insect SM, Ica f^m ■ 

Kay (1987) Um 

CAl pyramidal fj^a ■ 

Sala (1991) ttm 
SYMP N, 4mM Ca^+ /?„ ■ 

McCormick (1992) ttm 

TR N /3„ : 
Hutchcon (1994) 
TR N med dors 
Rybak (1997) 
Respiratory N 



0.013(y+20)~ 

■ i_e-(V'+20)/3 

0.31e-(^+20)/25 

1-6 

l_e-0072(l'-5) 

0.02(^+8.69) 

e(V+8.69)/5.36_]^ 

0.058(11. 3-y) 

■ e(ll-3-V)/i3.7„x 
0.085(V"+15.4) 

e{V+15.4)/9.9_l 
1.6 

0.02(y-1.31) 

e{V-1.3)/5.36)_l 



am+/3rj 



1 



>^m-ri-^m 



0.33 

Otrn+Pn 



Lindblad (1996) 



Rabbit AM 



Athanasiades (2000) 
MR 



Rhodes (2001) 
CORT pyramid 
Poirazi (2003) 
CAl pyramidal, soma 
Kager (2007) 
CA3 pyramidal 
Pospischil (2008) 
TR, CORT 



Pm ■ 

ah -- 

Ph-- 

a 

B* : 
* 

/3^ 

a-m 

Pm ■ 

Oim 

Pm 

O-m 

Pm 

a-m 

Pm 

ah -- 

Ph-- 



. 16.720(^+35.0) . 

■ l_e-(V+35.0)/2.5 + 

50 



-V/4.808 

_ 4.480(y-5.0) 

" e(^-5-0)/2-500_]^ 

8.490(V+28.0) 

g(V+28.0)/4.000_]^ 

67.922 

X_|^e^(^+28-0)/4.000 

= 0.346e°-°^25(y+20) 
= 1.8818e^o.oo732(y+20) 
0.2419e0-i45(y+20) 
0.0434e-°-°20i3(y+20) 

- 1600 

" i+e-(^-5)/l3-9 
20(V+8.9) 



- e(^+8-9)/5)-l 
_ 0.055(y+27.01) 

" l_e-Cl'+27.01)/3.8 

: 0.94e-(^+^3-°i)/^^ 

5.23(\/-71.5) 

" l_e-(V+''''71.5)/10 

0.055(V+27) 

" l_e-(V'+27)/3.8 

: 0.94e-(^+75)/i7 
0.000457e-(^+i3)/50 

0.0065 

l+e-e^+is)/^^ 



0.1 + 



1 



0.2 



am+/3n 



a + be ^ d > 



a=1.5, b= 211 

c=-37.427 

d=20.213 



17 + 



1 



<+Pl 



1 



ah+Ph 
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Table A3 - Experimental 
Time constants for L-type Ca^+ activation and, if included, VDI 

All cell types 



Source, cell type 


'Tm,max 


T/ 

''Tm,max 


Th,max 




T 7- 


Standcn (1982) 


19.82 


-28 








Insect SM, I^a 












Kay (1987) 


1.5-2.07 


-18 








Insect SM, Ipa 












Belluzzi (1991) 


0.564 


-8.0 


29.90 




-22 


SYMPN 












Sala (1991) 


1.16 


1 








SYMPN 












Thompson (1991) 






100 






CAl pyramidal 












Avery (1996) 


2.3 










CAS pyramidal 












Molitor (1999) 


1.5 


-10 








DCN N 












Donohoe (2000) 






Tfast = 17, 


Tslow = 57 


at mV 


Rat VM 






Tfast = 15, 


Tslow = 46 


RF 


Muinuddin (2000) 


4.0 ms distal 




n = 4402, 


T2 = 255 




Cat SMM, 20mM Ba^^ 


4.2 ms proximal 




Ti = 3905, 


T2 = 222 




Schnee (2003) 


« 0.5 










Turtle AHC 












Jackson (2004) 


« 2 










Dorsomcdial SCN 












Luin (2008) 


young 53.6 










Cultured human SM 


old 62.5* 
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Table A3 - Modeling 
Time constants for L-type Ca^+ activation and, if included, VDI 

All cell types 



Source, cell type 


Tm,max 


^Tm,max 


Th,max 


I' T^ ,max 


Traub (1991, 1994) 


2.11 


-18 






CAS pyramid 










McCormick (1992) 


1.541 


-15 






TR neuron 










Luo (1994) 


2.23 


-10 


25.2 


-10 


VM 






Min not max 




Hutcheon (1994) 


1.541 


-15 






TR med dors 










Jaffc (1994) 


1.5 








CAS pyramidal 










Wallenstcin (1994) 


1.541 


-15 






NRTN 










Migliore (1995) 


4.56 


-11 






CAS pyramidal 










Li (1996) 


0.134 


4 






DA SN 










Lindblad (1996) 


4.7 


-39 


226 


-37.4 


Rabbit AM 










LeBcau (1997) 


27 (all V) 








PC 










Booth (1997) 


40 (all V) 








Spinal MN 










Amini (1999) 


19.5 


-45 


GDI only 




MID DA 










Athanasiades (2000) 


0.563 


-28 


29.9 


-42 


MR 










Zhang (2000) 


3.8 


-37 


40 


-28 


Rabbit SA PM 






~ const, y > -28 


Shorten (2000) 


11.32 


-12 






PC 










Carlin (2000) (HVA) 


20, all V 


-10 






Carlin (2000) (LVA) 


20, all V 


-30 






Spinal MN 










Rhodes (2001) 


2.1 


-18 


20.0 




CORT pyramid 










Fox (2002) 


19.3 


-33 


230 


«< -60 


VM 
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Table A3 - Modeling, continued 



Source, cell type 


Tm,max 


^Tm,max 


Th,max 


^Tfi,max 


Poirazi (2003) 


0.826 


-35 


4.8 




CAl pyramidal soma 










Komendantov (2004) 


19.5 


-45 


GDI only 




MID DA 










Shannon (2004) 


2.38 


-15 


25.2 


-15 


VM 






Min not max 




Rhodes (2005) 


2.1 


-19 


200, all V 




TRN 










Bui (2006) 


20, all V 








Spinal MN 










Komendantov (2007) 


0.9, all V 








Mag Endocrine 










Pospischil (2008) 


7.04 


-38 


449.8 


-67 


TR, CORT 










Kapela (2008) 


3.65 


-40 


110 


-35 


Mesenteric SMM 
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Table A4 
Time constants if given by explicit formulas 

All cell types 



Source, cell type t„ 



Th 



Belluzzi (1991) 

SYMP N, 5mM Ca^+ 
Luo (1994)"^ 

VM 

Li (1996) 

DA SN 

Amini (1999) 
Komcndantov (2004) 
MID DA 
Shorten (2000), PC 

Fox (2002), VM 

Shannon (2004) 

VM 

Rhodes (2005) 

TRN 

Kapela (2008) 
Mesenteric SMM 



0.1+ ( 0.343e°-^25y *+i.88e 



-0.00732V 



l_e-(y+l0)/6-4 

0.035(y+10)(l+e-("^'+l0)/6'4) 

0.4 

e = -(y + ii)/8.3 



1.5 + 18e^ 



-{m^Y 



27 

g(V+60)/22_|_2e-2(V+60)/22 

0.25e-0-0iy o.o7e-0.05(V+40) 
IJ^e^O.mV "T ]^_|_g0.05(V+40) 

0.035(y+14.5)(l+e-(^+i4-5)/6) 



0.6 + 



cosh[(V+19)/24]) 



I V+40 ^ 

1.15 + e V 30 ; 



a+6e-(0-0337(V+lO))2 

a = 0.02, 6 = 0.0197 



30 + 



200 

l+e(v+20)/9.5 

1 



a+be-(0-0337(v+i4.5))^ 
a=0.02, b=0.0197 
200 



(^ y+35 ' 
45 + 65e ^ 25 > 
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Table A5 

Magnitude of L-type Ca^"*" currents or related quantities 

All cell types 



Source, cell type 


L-type magnitude 


Remarks 


Penington (1991) 


~ 4% of total Ica 


Soma, dissociated cells 


DRN SE 






McCormick (1992) 


- IpA/pF *, 


Dissociated cells 


TRN 


PL= 80.10-6 cm^/sec 




Schild (1993) 


QL = 1.44 mS/cm^ 


Cell, GL=0.075/iS, 


Rat med 




C=0.052 nF 


Hutcheon (1994) 


1.3mS/cm2 


PL as in McCormick(1992) 


TR med dors 






Kuryshev (1995) 


52% of HVA is L-type 


-70pA 


Rat corticotropes 




lOmM Ba?+ 


Li (1996) 


0.19 mS/ cm2 




DASN 


(dendrites) 




Booth (1997) 


0.33 mS / cm2 




Spinal MN 


(dendrites) 




Schroder (1998) 


13.2 fA (control) 


Ensemble averages 


Human VM 


38.2 fA (failing heart) 


single channel (Ba^+) 


Amini (1999) 


Gi=0.0005^S 


~ 30pA; gL=0.018 mS/ cm^ 


MID DA 




(sphere, rad 15/U ) 


Van Wagoner (1999) 


9.13 pA/pF 


3.35 pA/pF 


Human AM 


(control) 


(atrial fibrillators) 


Athanasiadcs (2000) 


Gi=0.075^S 


gL=1.44 mS/cm2 


MR 




cf Schild (1993) 


Carhn (2000) 






Spinal MN 






LVA & HVA, soma 


lOmS/ cm2 




LVA & HVA, dend 


0.3mS/ cm2 




Donohoe (2000) 


16.8 pA/pF 


at -10 mV 


Rat VM 


18.3 pA/pF 


RF 


Joux (2001) 


19.3 pA/pF 




SON 






Fox (2002) 


-- IpA/pF, 




VM 


P(^^ = 0.000026 cm/ms 




Mangoni (2003) 


3.32 pA/pF 


1.05 pA/pF 


Mouse SA PM 


(wild type) 


(1.3 knockout) 


Poirazi (2003) 


7 mS / cm2 




CAl pyramidal 


(soma) 
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Table A5 continued 



Source, cell type L-type magnitude Remarks 



Tipparaju (2004) 
Human AM 
Durante (2004) 
Rat DA SNc 
Komcndantov (2004) 
MID DA 
Rhodes (2005) 
TR N soma 
Vignali (2006) 
Pancreatic A,B 
Komendantov (2007) 
Mag Endocrine 
Marcantoni (2007) 
Adrenal CH cells 
Kapela (2008) 
Mesenteric SMM 
Xiang (2008) 
CAl pyramidal 
PuUer-Bicer (2009) 
Mouse VM 
Benitah (2010) 
VM 

Empson (2010) 
Purkinje neuron 



2.5 pA/pF 

Young adults 

27% of Ica is L-type 

0.216 mS / cm^ 

4 mS / cm^ 

60% of Ica is L-type 

0.4 mS / cm^ 
(soma) 

10% of Ica 

-60 to -40 mV 

100% of Ica is IcaL 

18.7 pA/pF 
(control) 
5.76 pA/pF 

WT 

Ica is mainly IcaL 

< 5% of Ica 



1.2 pA/pF, 2.6 pA/pF 
Infants, Older Adults 



0.22, 0.46 mS / cm^ 
(prim., sec. dendrites) 
~ 15 pA 



Ica is mainly IcaL 



3.17 pA/pF 
a2/S-l KO 



Soma &: dendrites 
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Table A6 
Parameters for examples of L-type channel subtype activation and VDI 



Source 


Subtype 


V^ 1 , Carrier 

'"■I- 


kfn 


^..^ 


kh 


Koschak (2001) 


Ca^l.2 


-3.9 


7.8 


-27.6 


13.8 


tsA201 




15 mM Ba2+ 










Cavl.3 


-17.5 

15 mM Ba2+ 


8.9 


-42.7 


6.6 


Xu (2001) 


Cavl.2 ai 


3.3 








Xenopus oocytes 


Cavl.2 ai 
Cavl.3 ai 
Cavl.3 ai 


5 mM Ca2+ 

-8.8 

5 mM Ba2+ 

-20 

5 mM Ca2+ 

-36.8 

5 mM Ba2+ 








Schnee (2003) 


Cavl.3 


-35 


4.7 


high frequ 




Turtle AHC 




2.8 mM Ca2+ 










Cavl.3 


-43 

2.8 mM Ca?+ 


4.2 


low frequ 




Lipscombc (2004) 


Cavl.2 


-16.1 


6* 






tsA201 


Cavl.3 


-40.4 

2 mM Ca2+ 


5.5* 






Hildebrand (2004) 


Cavl.2 


-13.9 




-33.7 


control 


tsA201, 2 mM Ba^+ 








-52.5 


allethrin 


Helton (2005) 


Cavl.2 


-17.6 








tsA201 


Cavl.3 


-39.4 

2 mM Ca2+ 








Cattcrall (2005) 


Cavl.l 


8 to 14 

10 mM Ca2+ 




-8 






Cavl.2 


-17 




-60 to -50 




HEK 




2mM Ca2+ 




-60 to -50 






Cavl.3 


-37 




-43 to -36 




HEK 




2 mM Ca?+ 




-43 to -36 






Cavl.4 


-2.5 to -12 




-27 to -9 




HEK 




2 - 20 mM Ca2+ 




-27 to -9 




Navcdo (2007) 


Cavl.2 


31.9 

110 mM Ba?+ 


9.3 


-2.7 


5.2 


tsA201 


Cavl.3 


9.7 


10.8 


-18.7 


3.5 


Arterial SMM 




33.8 


8.0 


-4.9 


5.4 


Shapiro (2007) 


Cavl.2 


-20 


6 






Model, MN 


Cavl.3 


-41 


6 
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Table A6 continued 



Source 


Subtype 




V^ 1 , Carrier 

'"■'9. 




kfn 


^M 


kh 


Luin (2008) 


Ca^l.l 




8.73 young 




7.72 


-5.77 


7.30 


Cultured human SM 


(Assumed) 




10 mM Ca2+ 












Ca^l.l 




14.47 old 




7.74 


2.38 


10.1 


Peloquin (2008) 
tsA201 


Cavl.4 




5 

20 mM Ba2+ 

-13 

20 mM Ba2+ 


23°C 

37° C 


5.9 
3.9 






Singh (2008) 


Cavl.3 (42, long) 
Cavl.3 (42A, short) 

Cavl.3 (42, long) 
Cavl.3 (42 A, short) 


-2.2 

-12.9 

15 mM Ca2+ 

-11.8 

-20.2 

15 mM Ba2+ 




9.1 
6.9 

7.8 
6.6 






FuUer-Biccr (2009) 


WT 




-3.5 




4.87 


-15.15 


5.54 


Mouse VM 


02/(5-1 KO 




4.66 

1.8 mM Ca2+ 




5.88 


-2.18 


6.14 


Putzier (2009a) 


Cavl.3 




-31.1 




5.35 






Cordeiro (2009) 


Cavl.2 /32 WT 


1.5 






-24.8 


9.13 


tsA201 


B mutant 




3.5 

15 mM Ca2+ 






-30.0 


6.25 


Gebhart (2010) 
tsA201 


Cavl.3 /33 




-0.9 

15mM Ca2+ 




9 
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